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"N  A  TMq'  mode  is  incident  in  a  circular  waveguide  in  which  only  the  TE^  and 

TM_.  modes  can  propagate.  Two  apertures  in  this  circular  waveguide  feed  two 
rectangular  waveguides  of  identical  cross  section  in  which  only  the  TE- _  domi¬ 
nant  mode  can  propagate.  ^^le  close  each  aperture  with  perfect  condictors. 
Attempting  to  restore  the  original  tangential  electric  field  on  both  sides  of 
each  of  these  conductors,  wp  pla™^  the  magnetic  current  M  ,on  one  side  of  the 

conductor  and  -M  on  the  other  side  of  the  conductor.  We  expand  M  as  a  linear 
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20.  ABSTRACT  (continued) 

combination  of  expansion  functions  each  of  which  is  approximately  the  magnetic 
current  associated  with  one  of  the  rectangular  waveguide  modes.  Substituting 
this  ^into  the  equation  of  continuity  of  the  tangential  magnetic  field  across 
the  apertures  and  testing,  as  in  Galerkin’s  method,  this  equation  with  the 
expansion  functions  for  I£,  -w^arrive^at  a  matrix  equation  which  determines  the 
coefficients  of  the  expansion  functions  in  the  expansion  for  M.  In  turn,  M 
determines  the  field  in  the  waveguides,  and,  in  particular,  the  amplitudes  of 
the  TE-^q  modes  in  the  rectangular  waveguides.  Detailed  formulas  for  the 
matrix  elements  are  given  in  terms  of  Bessel  functions  and  roots  of  Bessel 
functions.  Further  work  is  required  to  program  these  formulas  in  order  to 
obtain  numerical  results  using  a  digital  computer. 
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Chapter  1 

Statement  of  the  Problem 


Two  apertures  in  a  circular  waveguide  feed  two  identical  rectangular  wave¬ 
guides  as  shown  in  Fig.  1.1.  The  walls  of  each  waveguide  are  perfectly  con¬ 
ducting.  The  interiors  of  the  left-hand  waveguide,  the  right-hand  waveguide, 
and  the  circular  waveguide  axe  called  regions  1,  2,  and  3,  respectively.  Ho¬ 
mogeneous  space  of  permeability  p  and  permittivity  e  exists  in  regions  1,  2, 
and  3.  The  excitation  is  a  TMo\  wave  of  unit  amplitude  traveling  in  the  z 
direction  in  the  circular  waveguide.  The  circular  waveguide  is  of  radius  a  and 
is  terminated  by  a  perfectly  conducting  wall  at  z  =  L3.  The  radius  a  is  such 
that  only  the  TE\\  and  T  Mq\  modes  can  propagate  in  the  circular  waveguide. 
A  problem  similar  to  the  one  being  described  was  previously  treated  in  [1]. 

Both  rectangular  waveguides  run  parallel  to  the  x  axis.  Both  have  the 
same  cross  section  (— |<y<|,— |<z<|)  where  c  <  b  and  b  is  such  that 
only  the  TE\0  dominant  mode  can  propagate  in  each  rectangular  waveguide. 
The  aperture  which  feeds  the  left-hand  rectangular  waveguide  in  Fig.  ]  1  is 
called  A\.  This  aperture  is  the  surface  for  which  (p  =  a,  -k  —  <p0  <  <t>  < 
7r  4-  0O,  -f  <  z  <  §)  where 

p  =  \]x2  +  y2  (1.1) 

<t>o  =  sin_1(^)  (1.2) 

The  aperture  which  feeds  the  right-hand  rectangul?:  waveguide  in  Fig.  1.1 
is  called  A?.  This  aperture  is  the  surface  for  which  (p  =  a,  —  <p0  <  (p  < 

-f  <  *  <  f). 
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Figure  l.l:  Top  and  side  views  of  the  T A/0i  to  T El0  mode  converter. 
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The  voltage  to  current  ratio  of  the  TE\q  mode  in  region  1  is  taken  to 
be  Z\  at  x  =  —  L\.  All  other  rectangular  waveguide  modes  are  evanescent. 
The  voltage  to  current  ratios  of  the  evanescent  modes  at  x  =  —  Lx  do  not 
come  into  play  because  L\  is  taken  to  be  so  large  that  any  evanescent  wave 
emanating  from  the  termination  at  x  =  —  Lx  will  have  negligible  amplitude 
upon  arrival  at  the  apertur®  A\.  The  voltage  to  current  ratio  of  the  T Ew 
mode  in  region  2  is  taken  to  be  Z2  at  x  =  L2.  Here,  L2  is  taken  to  be  so 
large  that  any  evanescent  wave  emanating  from  the  termination  at  x  =  L2 
will  have  negligible  amplitude  upon  arrival  at  the  aperture  A2. 

As  previously  stated,  the  excitation  is  a  z  traveling  (traveling  in  the  2 
direction)  TM01  wave  of  unit  amplitude  in  the  circular  waveguide.  By  as¬ 
sumption,  this  is  the  only  z  traveling  wave  at  z  =  —  |.  This  wave  is  produced 
by  an  impressed  source  or  sources  located  in  the  region  for  which  z  <  za  <  —  | 
of  the  circular  waveguide.  Since  neither  of  the  rectangular  waveguides  nor 
the  region  for  which  z  >  za  of  the  circular  waveguide  contains  an  indepen¬ 
dent  source,  the  electromagnetic  field  for  z  >  z,  in  the  circular  waveguide 
depends  only  on  the  amplitudes  of  the  z  traveling  waves  at  z  =  z,  and  not 
on  the  manner  in  which  these  amplitudes  were  produced.  Therefore,  we  can, 
without  loss  of  generality,  assume  that  the  impressed  source  of  the  unit  am¬ 
plitude  z  traveling  TM0 1  wave  in  the  circular  waveguide  is  an  electric  current 
source  ^imP  whose  —  z  traveling  (traveling  in  the  —z  direction)  waves  see  a 
matched  load,  that  is,  these  waves  are  never  reflected.  As  shown  in  Fig.  1.1, 
limp  is  located  at  z  <  —  |  in  the  circular  waveguide.  The  objective  is  to  find 
the  electromagnetic  field  in  regions  1  and  2  of  Fig.  1.1  and  in  the  portion  of 
region  3  for  which  z  >  z,  in  Fig.  1.1. 
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Chapter  2 
Formulation 


Following  the  generalized  network  formulation  for  aperture  problems  [2],  [3], 
we  close  the  apertures  A\  and  A2  with  perfect  electric  conductors  of  infinites¬ 
imal  thickness.  As  shown  in  Fig.  2.1,  we  place  the  surface  density  of  magnetic 
current  on  the  region  1  side  of  the  closed  aperture  Ax,  — on  the 
region  3  side  of  A\,  on  the  region  2  side  of  the  closed  aperture  A2,  and 
— M ^  on  the  region  3  side  of  A2.  In  Fig.  1.1,  the  tangential  electric  field 
and  the  tangential  magnetic  field  are  continuous  across  Ax  and  A2.  The 
arrangement  of  magnetic  currents  in  Fig.  2.1  ensures  continuity  of  the  tan¬ 
gential  electric  field  across  A\  and  A2.  Now,  the  fields  in  Fig.  2.1  will  be  the 
same  as  those  in  Fig.  1.1  if  M .  and  M are  adjusted  so  that  the  tangential 
magnetic  field  in  Fig.  2.1  is  continuous  across  A\  and  A2. 

Continuity  of  the  tangential  magnetic  field  across  Ax  in  Fig.  2.1  is  ex¬ 
pressed  as 

(2.1) 

where  is  the  magnetic  field  in  region  1  and  is  the  magnetic  field 
in  region  3.  In  (2.1),  the  subscript  “tan”  denotes  the  components  tangent 
to  A\.  Continuity  of  the  tangential  magnetic  field  across  A2  in  Fig.  2.1  is 
expressed  as 

*2o  =  4a„  »■>  ^  (2-2) 

where  is  the  magnetic  field  in  region  2.  In  (2.2),  the  subscript  “tan” 
denotes  the  components  tangent  to  A2. 

The  electromagnetic  field  in  region  1  of  Fig.  2.1  is  due  to 
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Figure  2.1:  Top  and  side  views  of  the  situation  equivalent  to  that  of  Fig.  1.1 
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(2.3) 

(2.4) 

The  superscript  “(1)”  on  R  and  R  on  the  right-hand  sides  of  (2.3)  and  (2.4) 
denotes  radiation  in  region  1  of  Fig.  2.1.  The  UQ”  on  the  right-hand  sides  of 
(2.3)  and  (2.4)  indicates  that  there  is  no  electric  current  source. 

The  field  (£*2\i£^)  in  region  2  of  Fig.  2.1  is  due  to  M ^  placed  on  the 
region  2  side  of  the  closed  aperture  Ay. 

R(2)  =  £W(^M.{2])  (2.5) 

R(2)  =  £(2>(Q,  M.{2))  (2.6) 

The  superscript  “(2)”  on  R  and  K  on  the  right-hand  sides  of  (2.5)  and  (2.6) 
denotes  radiation  in  region  2  of  Fig.  2.1.  The  “Q”  on  the  right-hand  sides  of 

(2.5)  and  (2.6)  indicates  that  there  is  no  electric  current  source. 

The  field  (R^3\H_^)  for  z  >  za  in  region  3  of  Fig.  2.1  is  due  to  the  current 
sources  J}m P,  —M}1\  and  —  radiating  in  the  circular  waveguide  with 
the  apertures  Ai  and  A<i  closed  by  perfect  conductors,  with  the  perfectly 
conducting  wall  at  z  —  L$,  and  with  a  matched  load  at  z  =  zl  where  zL  is 
any  value  such  that  all  of  j}mP  lies  in  the  region  for  which  z  >  zi\ 

£(3)  =  £(3)(  J'imp^  _  £<3)(Q?  ^£(1))  _  £(3)(Q)  (2.7) 

£(3)  =  £(3)(2im P  q)  _  £(3)(Q?iJ£ri))  -£(3)(Q,M(2))  (2.8) 

The  superscript  “(3)”  on  the  right-hand  sides  of  (2.7)  and  (2.8)  denotes 
radiation  in  region  3  with  the  apertures  closed,  with  the  short  at  2  =  L3, 
and  with  the  matched  load  at  the  other  end.  The  first  argument  of  each  field 
on  the  right-hand  sides  of  (2.7)  and  (2.8)  is  an  electric  current  source;  the 
second  argument  is  a  magnetic  current  source. 

Substitution  of  (2.4)  and  (2.8)  into  (2.1)  gives 

-rC(Q,  At"1)  -flgjll,  M"»  )-iC(2.  M.m) 

Substitution  of  (2.6)  and  (2.8)  into  (2.2)  yields 

-ttgn(Q,M",)-aS„(a,M<I,)-itSn(Q.M<!|) 

(2.10) 


= -itgn UimP,Q)  on /l, 

(2.9) 

= -Cn«imp'2> 


placed  on  the  region  1  side  of  the  closed  aperture  Ay. 

£ri)  _  £(1)(Q,M(1)) 

2£(1)  =  £(1)(Q,M(1)) 
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pTM  QTM  pTE  qTE 

m(1)  =EE  v™m£m(4>,z)  +  £  E  vr%T(^>  (2.11) 

p=l  7=1  p=0  7=0 

P+9^0 

pTM  QTM  pTE  qTE 

m®  =  E  E  (£r"JC'<#,*)  +  E  E  (2.12) 

p=l  7=1  p=0  7=0 

P+99^0 

where  the  V’s  axe  unknown  coefficients  to  be  determined  and 
Mj*(<M)  =  J4efM(y1+,x+)  +u,S-~^esypq(y1+,z+), 


P  —  a 

n-<l>o<<f><'x-lr4>o 
_£  <  z  <  £ 

2  —  2  —  2 


(2.13) 


40<M)  =  34efp<J(y2+,x+)  -iit^^eJM(t/2+,2+),  j  -<i>0  <  (j>  <  <t>a  (2.14) 

<p0  I  _£  <  z  <  £ 

2  —  —  2 


o  =  a 


where  and  u*  are  the  unit  vectors  in  the  <f>  and  z  directions,  respectively. 
Furthermore,  and  e6^  are,  respectively,  the  y  and  2  components  of  £%q. 
Now,  6  is  either  TM  or  TE,  and  §££*  and  e™  are  given  by  (A.  10)  and  (A. 23), 
respectively.  In  (2.13)  and  (2.14), 


y1+  =  (*•-  #c<>+  | 
y2+  =  d>x0  +  | 
z+  =  2  +  | 


(2.15) 

(2.16) 
(2.17) 


where 


a  sin  4>0 


(2.18) 


In  Chapter  3,  it  is  shown  that  (ff(1)(Q,  ft/p™).  ff(1)( 0.  Afp™))  is  approxi¬ 
mately  a  TMp,  field  in  region  1,  that  (.F^fO.  ff (1)(0.  MlEE))  is  ap¬ 

proximately  a  TEpq  field  in  region  1,  that 

is  approximately  a  TMpq  field  in  region  2,  and  that  ( E^H 0.  M^E),  H^iO, 
M%E))  is  approximately  a  T  Epq  field  in  region  2. 


Substituting  (2.11)  and  (2.12)  into  (2.9)  and  (2.10),  we  obtain 

pTM  qTM 

-  E  E  C"UC»-^™)  +  £tan(^™)} 

P=l  <7=1 
pTE  qTE 

-EE  l(a,JO+flgn(&*0) 

7=0  <7=0 

p+9^0 

pTM  qTM  pTE  qTE 

-E  E 

P=1  <7=1  p=0  9=0 

P+9/0 

=  -flgnUin>P,S)  on  A,  (2.19) 


pTM  qTM  pTE  qTE 

-EE  iff-fOi.  JO  -  E  E 

P=1  9=1  p=0  g=0 

P+9*° 

pTM  qTM 

-EE  v^«(i£«!i1(Q,<-“)  +  i£gn(fl,M;™); 

P=1  9=1 

pTE  qTE 

-EE  C*<flSn<a..40+iC(2-J0) 

p=0  9=0 
P+9^0 

=  -i£^(iimp,!l)  on  /t2  (2.20) 

The  symmetric  product  <  A,  B  >  between  two  vectors  A  and  B_  is  defined 
to  be  the  surface  integral  of  their  dot  product  over  whichever  aperture  they 
axe  defined: 

<  A,B  >=  [  f  A  Bdis  (2.21) 

J  Jai  or  a2 

Here,  ds  is  the  differential  element  of  surface  area.  First  taking  the  symmetric 
product  of  (2.19)  with  each  of  the  expansion  functions  M ^Tn'w  and  M.^n  that 
appear  in  (2.11),  and  then  taking  the  symmetric  product  of  (2.20)  with  each 
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of  the  expansion  functions  M 2™  and  M2™  that  appear  in  (2.12),  we  obtain 


[y1  +  y2  +  y3] 


■  ylTM  ■ 

’  f\TM  ' 

ylTE 

flTE 

y2TM 

pTM 

y2TE 

pTE 

(2.22) 


— • 

where  the  V’s  are  column  vectors  of  the  unknown  V’s  in  (2.11)  and  (2.12). 
— • 

The  P s  axe  also  column  vectors;  the  right-hand  side  of  (2.22)  is  called  the 
excitation  vector. 

The  elements  of  the  V’s  in  (2.22)  are  given  by 


_  yi& 

vpq  ) 


j  =  1,2,  •  •  •  iV5 
<  7  =  1,2 
8  =  TM,  TE 


(2.23) 


The  subscript  “j”  on  the  left-hand  side  of  (2.23)  is  a  condensation  of  the 
double  subscript  “jxj”  on  the  right-hand  side  of  (2.23);  a  one  to  one  corre¬ 
spondence  is  established  betweeen  each  pair  of  integers  p  and  q  in  use  and 
each  of  the  positive  integers  1,2,  •  ■  ■  N6 .  The  correspondence  in  the  “TE” 
expressions  (8  =  TE  in  (2.23))  may  be  different  from  that  in  the  “TM” 
expressions  (6  =  TM  in  (2.23)).  The  elements  of  the  7’s  in  (2.22)  are  given 
by 


If  =  -  j jA  ■  i£(3)Uimp,  Q)  <«», 


'  i  =  1,2,  •  •  • , 

<  Q  =  1,2 

0  =  TAf,  TE 


(2.24) 


where  i  is  related  to  mn  in  the  same  way  as  j  is  related  to  pq  in  (2.23).  The 
relationship  in  (2.23)  was  of  the  type  “(5”;  that  in  (2.24)  is  of  the  type  "3'' . 
The  subscript  “tan”  which  was  attached  to  Q)  in  (2.19)  is  not 

needed  in  (2.24)  because  M°gn  is  tangent  to  Aa. 

In  (2.22),  y1,  y2,  and  Y 3  are  the  admittance  matrices  for  regions  1,  2. 
and  3,  respectively.  The  matrix  y1  is  given  by 


y1 


yl,lTM,lTM  yl,lTM,lTE 
yl,\TE,lTM  y\,\TE,\TE 

0  0 


0 

0 

0 


0  ' 
0 
0 
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(2.25) 


where  each  entry  of  the  axray  on  the  right-hand  side  is  a  submatrix.  Each 
element  of  the  “0”  submatrices  is  zero.  The  elements  of  the  non-zero  subma¬ 
trices  are  given  by 

YiMU  =  - ■  j0==TMlTEN‘  (2'26) 

8  =  TM,TE 

where  j  is  related  to  pq  as  in  (2.23);  i  is  similarly  related  to  mn.  Nevertheless, 
the  relationship  between  i  and  mn  may  be  different  from  that  between  j  and 
pq.  The  relationship  between  i  and  mn  is  of  the  type  The  relationship 
between  j  and  pq  is  of  the  type  US”.  The  matrix  Y 2  is  given  by 

"0  0  0  O' 

y2  _  ®  0  0  0  .  9 n\ 

1  —  Q  Q  Y2’2™’2™  Y2'2TM,2TE 

Q  Q  Y2’2TE'2™  Y2'2TE'2TE 


Y2  = 


(2.27) 


where  the  elements  of  the  “0”  submatrices  are  zero.  The  elements  of  the 
non-zero  submatrices  are  given  by 


=  Km(  a,  m%)  is. 


V« 

j  =  1,2,  •  •  •  ,iV5 
f3  =  TM,TE 
8  =  TM,  TE 


(2.28) 


where  i  is  related  to  mn  and  j  is  related  to  pq  as  in  (2.26).  The  matrix  Y 3 
is  given  by 

-  y3’1™’1™  y3’iTM'iTE  y3,1™'2™  y3,1™'2TE  ' 

y3,lXE,lTAf  y3,lTE,lTE  y3,lTE,2TM  y3,irE,2TE 

Y  —  y3'2™'1™  y3'2™'iTE  y3'2™'2™  y3'2tm,2TE  (2.29) 

y3,2TE,\TM  y  3,2TE,lTE  y3,2TE,2TM  y3.2TE.2TE 


(2.29) 


where 


r3,a0,yS 


=  -JJAM°l-Km(a,M£)is, 


j  =  1,2 
a  =  1,2 

0  =  TE,  TM 

7  =  1/2 

8  =  TE,  TM 


(2.30) 
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In  (2.30),  i  is  related  to  mn  and  j  is  related  to  pq  as  in  (2.26). 

If  the  elements  of  the  admittance  matrices  Yl,  Y2,  and  Y3  and  the  column 
vectors  P™,  PTE,  f1™ ,  and  PTE  can  be  evaluated,  then  (2.22)  can  be 
solved  for  V1™,  VlTE,  V2™,  and  V2TE.  These  V’s  determine  and 
according  to  (2.11)  and  (2.12).  Next,  M and  can  be  substituted 
into  expressions  (2.3)-(2.8)  for  the  fields  in  the  waveguide  regions. 
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Chapter  3 

The  Admittance  Matrices  for 
the  Rectangular  Waveguides 


In  this  chapter,  Y>f0'xS  of  (2.26)  and  Y2-20'26  of  (2.28)  are  evaluated.  We 
approximate  Y*'10,1S  of  (2.26)  by  Yt]’10'16  given  by 


yl,10,lS 


i  =  l,2,---,N0 
j  =  1»2,  ■  •  •  ,7V5 
0  =  TM,TE 
.  8  =  TM,TE 


We  approximate  Y2-'20,26  of  (2.28)  by  Y2'20'26  given  by 


y  2, 2/3, 25 

rij 


J m!1  *  K{2){0.,  M.2pq)ds, 


j  =  1, 2,  -  •  • ,  iV6 
0  =  TM,TE 
6  =  TM ,  TE 


(3.1) 


(3.2) 


In  (3.1)  and  (3.2),  M™  and  M™  are  approximations  to  M**  of  (2.13)  and 

a  25  *  * 

Mpq  of  (2.14).  Furthermore,  A\  and  A2  are,  respectively,  the  surfaces  on 

*  15  *25 

which  Mvq  and  Mvq  are  defined.  In  (3.1)  and  (3.2),  j  is  related  to  pq  as  in 
(2.23);  i  is  similarly  related  to  mn. 

The  Y' s  of  (3.1)  and  (3.2)  will  be  evaluated  by  first  defining  the  M’s  that 
appear  in  (3.1)  and  (3.2)  and  then  by  finding  the  i£’s  that  appear  in  (3.1) 
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and  (3.2).  The  electric  fields  associated  with  these  2£’s  will  also  be  found 
because  these  electric  fields  are  needed  in  order  to  evaluate,  as  indicated  in 
the  last  two  sentences  of  Chapter  2,  the  field  in  the  rectangular  waveguides. 
Finally,  the  M’s  and  the  i£’s  will  be  substituted  into  (3.1 )  and  (3.2)  to  obtain 
appropriate  expressions  for  the  Y's. 

We  define  by 


=  M^q{y,z)  =  -4<,(y+'z+)  x  M*  .  { 


x  =  -x0 

-§<*<! 


(3.3) 


where  z+  and  x0  are  given  by  (2.17)  and  (2.18),  respectively.  Furthermore, 


y+  =  V  +  f  (3-4) 

As  defined  by  (3.3),  to*  exists  on  a  portion  of  the  x  =  —  x0  plane.  With 

x0  given  by  (2.18),  —  xa  is  the  average  value  of  x  over  the  curved  surface  on 
. ,  » is 

which  Mpq  of  (2.13)  exists.  Mpq  is  placed  on  the  region  1  side  of  a  perfect 
conductor  which  covers  the  plane  surface  in  (3.3).  This  surface  is  called  Ax. 
We  relate  y  of  (3.4)  to  <j>  by 


y  =  (»  -  <t>)xa  (3.5) 

Relationship  (3.5)  is  reasonable  because: 

1)  y  =  0  when  <f>  =  x 

2)  y  =  |  when  <f>  =  -k  —  <j>0 

3)  y  =  —  2  when  <f>  =  7r  +  4>0 

The  above  items  2)  and  3)  are  obtained  by  using  (2.18)  and  (1.2).  Given 
(3.5),  the  interval  (y,  y  +  dy)  corresponds  to  the  interval  (</>,  <$>  +  d0)  where  o 
is  related  to  y  by  (3.5)  and 

d<p  =  (3.6) 

x0 

'  15 

Here,  dy  and  d<f>  are  the  differentials  of  y  and  <f>,  respectively.  M fq  of  (3.3)  is  a 
good  approximation  to  Mpsq  of  (2.13)  because  the  transverse  and  longitudinal 
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«  16 

voltages  produced  by  Mpq  on  the  portion  of  its  surface  for  y  in  the  interval 
(y  4-  dy)  and  z  in  the  interval  ( z  +  dz)  are  the  same  as  the  transverse  and 
longitudinal  voltages  produced  by  Mps;  on  the  corresponding  portion  of  its 
surface.  The  latter  portion  of  surface  is  the  surface  for  which  p  =  a,  <f>  is  in 
the  interval  ( <f> ,  <j>  +  d<f>),  and  z  is  in  the  interval  (z,  z  +  dz).  Here,  dz  in  the 
differential  of  z.  Moreover,  <f>  and  d(f>  are  related  to  y  and  dy  by  (3.5)  and 
(3.6),  respectively. 

To  show  equality  of  the  voltages  mentioned  in  the  previous  paragraph, 
»  15  . , 

we  note  that  Mpq  of  (3.3)  produces  the  tangential  electric  field  ej*  given  by 

%=-a (3.7) 

on  its  surface  and  that  Mpsq  of  (2.13)  produces  the  tangential  electric  field 
45  given  by 

45  =  x  =  -]^^~£e5P,(y1+,2+)  +  uzeszpq{yl+,z+)  (3.8) 

on  its  surface.  Here,  14,  uy,  u.z,  Si and  14  are  unit  vectors  in  the  directions 
indicated  by  the  subscripts.  The  transverse  voltage  produced  by  M  pq  is  the 
voltage  at  y  +  dy  with  respect  to  that  at  y.  This  voltage  is  Vt  given,  in  view 
of  (3.7),  by 

v  =  45  •  (-MsO  =  -elv^y+'z+)dy  (3-9) 

The  transverse  voltage  produced  by  Mlpq  is  the  voltage  at  (f>  +  d(f>  with  respect 
to  that  at  <j).  This  voltage  is  Vt  given  by 

Vt  =4 5  •  (-2 uQd<f>)  =  -t\pq  •  (y1+,  z+)dy  (3.10) 

Equations  (2.18),  (3.6)  and  (3.8)  were  used  to  obtain  (3.10).  The  longitudinal 

*  b  t 

voltage  produced  by  M^q  is  Vt  given  by 

Vi  =  45  •  (-Uzdz)  =  -ezpq(y+,z+)dz  (3.11) 

The  longitudinal  voltage  produced  by  Afpq  is  Vt  given  by 

Vt  =  4? '  (-M2)  =  -4pq(yl+^r)d2  (3.12) 
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From  (2.15),  (3.4),  and  (3.5),  we  have  y1+  =  y+  so  that  it  is  now  evident 

from  (3.9)— (3. 12)  that  the  transverse  and  longitudinal  voltages  produced  by 

•  is  -is 

Mpq  are  the  are  the  same  as  those  produced  by  Mpq. 

We  define  M™  by 

=  M^(y, 2 )  =  4q(y+,z+)  x  Ux  ,  J  -f  <  y  <  f  (3.13) 

l  -§<*<§ 

where  z+,  x0 ,  and  y+  are  given  by  (2.17),  (2.18),  and  (3.4),  respectively. 

of  (3.13)  is  a  good  approximation  to  Mps  of  (2.14)  because  the  transverse  and 

*  26 

longitudinal  voltages  produced  by  Mpq  on  a  differential  element  of  its  surface 
are  the  same  as  the  transverse  and  longitudinal  voltages  produced  by  Mpsq 
on  the  corresponding  differential  element  of  its  surface.  Equality  of  these 
voltages  follows  from  an  argument  similar  to  that  in  the  previous  paragraph. 

The  problem  of  finding  the  field  (£^(Q,  M**).  M**))  is  a  boundary 

value  problem  in  which  the  transverse  electric  field  at  x  =  —x0  in  region  1  is 
given  by  (3.7)  as 

£^  =  4,(5/\z+)  (3.14) 

Furthermore,  region  1  is  terminated,  as  in  Fig.  2.1,  at  x  =  —Li  with  the 
impedance  Z\.  Therefore,  as  explained  in  the  third  paragraph  of  Chapter  1, 
the  only  x  traveling  wave  at  z  =  —  x0  is  a  TE\$  wave. 

We  choose 

(£<l)(Q  (3.15) 

(£‘,,(Q,&!”),i£(,,(ll,M|”))  =  Cj0E-(Ejt,Rw-) 

+cr/+(£foE+.a.foE+)  (3.16) 

(£|I)(Q,jO,fl|,,(a,iO)  =  C™- (p,«  #  1.0)  (3.17) 

where  the  C’s  cure  unknown  constants.  The  (i£,  2£)’s  on  the  right-hand  sides  of 
(3. 15)— (3.17)  are  mode  fields  given  by  (A. 3),  (A. 14),  and  (A. 15).  Substituting 
(A. 3),  (A. 14),  and  (A. 15)  into  (3.15)-(3.17),  we  obtain 
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£(1>(Q,iO  =  c™~  [-ZjqM^M(y\z^)  +  ^ 


i£(1)(Q,C) 

kw&m£e) 


(3.18) 

C™~L™(y+,z+)S”x  (3.19) 

(C?E+e~™x  +  CXE~ enox)^Q  (y+ , z+)  (3.20) 

{C™+e~™x  -  C™-e™)Y™hZ£{y\z+) 
+U*(C?0E+e~™  +  —  (3.21) 


rw&m1te)  = 

i£(1)(Q,MiD  = 


Ctf-tfiy^z+W"1,  (?,?)#(  1,0) 

/~*TE—  f  \/TElTE/„.+  _+\  .  klqrfqE(y+,z+) 


C™~ 


P9  ^9 


(3.22) 


(3.23) 


Setting,  as  required  by  (3.14),  the  transverse  part  of  £(1)(Q,MP,)  equal 
to  fi£,(2/+,z+)  when  i  =  — 1„,  we  obtain 


rrM-  _  _£!!!: 

7TM 

P<7 


1  =  CjE+e™x°  +  Cl 'E~e~™x° 

C£E-  =  e’"‘-,  (P,?)^(1,0) 


(3.24) 

(3.25) 

(3.26) 


The  presence  of  Zx  requires  that,  at  x  =  —  Z,x,  the  ratio  of  the  coefficient  of 
£uf(y+>2+)  *n  (3.20)  to  the  coefficient  of  kJo{y+,z+)  in  (3.21)  be  — Zx: 


y  \rTE  _  ^10 


CTE+e-noLi  +CTf~e-™Ll 


(3.27) 


1  10  "  C?0E+e™^  -  Ctf-e-noLt  v^' ' 

The  constants  Cf0E+  and  CXE~  that  satisfy  (3.25)  and  (3.27)  are 

-  ^TE+  =  _ (Zx^-l)e-^ _ 

10  2  {sinh(7io(ii  -  x0))  +  ZXY^E  cosh(7x0(Ii  -  x0))} 


rTE- 
^  in 


_ {ZxYjE  +  l)e~”°Li _ 

2  {sinh(710(Ii  -  x0))  +  ZXY?QE  cosh(710(Li  -  x0))} 


(3.29) 


In  view  of  (A. 13)  and  (A. 25),  substitution  of  (3.24),  (3.26),  (3.28),  and 
(3.29)  into  (3. 18)— (3.23)  gives 


a">(a,  iC)  = 


7p<? 


7p,(x  +  x0) 


(3.31) 


£(1)(Q,Mi0  )  =  g0E(y+,z+) 

sinh  (7iq(Ii  +  x))  +  ZjY&E  cosh  (710 (L1  +  x)) 
sinh  (7io(Zi  -  x0 ))  +  Z\Y&E  cosh  (710(^1  -  x0 )) 


i£(1)(Q,Mio )  =  fi0E(y+,z+) 
i7io(cosb  (7iq(Ii  +  x))  +  ZxY?qe  sinh  (710(^1  +  x)))  _ 
u;/i(sinh  (710(Z,i  -  J0))  +  ZiY&E  cosh  (710(Li  -  x0)))  ~x' 

,te,  +  +,  J* i2o(sinh  (710(^1  +  x))  +  ZXY?0E  cosh  (710(^1  +  *)))  ^ 

10  ’  ’ufi{ sinh  (7io(A  -  *•))  +  ZiY&E  cosh  (710(^1  -  *0))) 


£(,)(a,  jO  =  47(»+ ,  *V”(l  +  x°] ,  (P, ,)  #  (1,  o) 


klqrPlE{y+,z+Y 


H.(X\<l,h£E)  =  ^  2&V\*+)  -H.-H2EI 

w/x  V  7p, 

7p,(x  +  x„) 

■«  1  0m)  #  l1-0) 


(3.34) 


(3.35) 


The  electromagnetic  field  in  (3.30)  and  (3.31)  is  a  —  x  traveling  wave  (A.3), 
that  in  (3.32)  and  (3.33)  is  a  combination  of  x  traveling  and  —  x  traveling 
waves  (A. 14)  and  (A. 15),  and  that  in  (3.34)  and  (3.35)  is  a  —  x  traveling 
wave  (A. 15).  When  x  =  —x0,  the  transverse  part  of  the  electric  field  (3.30) 
is  the  electric  field  (3.32)  is  g[E,  and  the  transverse  part  of  the  electric 
field  (3.34)  is  §£E.  When  x  =  —Lx,  the  ratio  of  the  TE\0  voltage  associated 
with  (3.32)  to  the  TE\0  current  associated  with  (3.33)  is  —  Zx. 

In  a  development  similar  to  that  in  the  previous  three  paragraphs,  we 
obtain 


.  a  2TM. 
£(2)( ) 


7p<J(x  j0) 
e 


(3.36) 
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(3.37) 


i£w(a ,jO  =  ^CV,*+)*  7”(i  *‘) 

7m 

£(3)(Q,MioE)  =  d?(»+.*+) 

sinh  (710(Za  ~  *))  +  cosh  (710(^2  ~  *))•  m 

sinh  (7io(Ia  -  *«))  +  Z2Y&E  cosh  (710(^2  -  *0)) 

R(2)(0.  ,M?oE)  =  -h™(y+,z+) 

J7io(cosh  (7io(Zr2  -  a))  +  Zrfjxjj sintt  (Wjg  ~  *)))  _ 
u;/i(sinh  (710(12  -  *«))  +  ^2^10^  cosh  (710(^1  -  x0)))  *** 

,T£/  +  +  x  jfe?0(sinh  (710(^2  -  a))  +  £2*10*  cosh  (710(^2  -  *)))  ,«  „Q) 

‘W 10  W  ’  V(sinh  (7io(£2  -  x.))  +  Z*Y&*  cosh  (710(^2  -  *•)))  1 ' 

£<>,(a,O-4*(»+.*+)‘"7"(:'"*').(p.«)>‘(i.0)  (3.40) 

Km%&CB)  =  (iffV,  *+)  +  M.  e_7w(l  '  I“), 

(p,g)/(l,0)  (3,41) 

The  electromagnetic  field  in  (3.36)  and  (3.37)  is  an  x  traveling  wave  (A.2), 
that  in  (3.38)  and  (3.39)  is  a  combination  of  x  traveling  and  — x  traveling 
waves  (A. 14)  and  (A. 15),  and  that  in  (3.40)  and  (3.41)  is  an  x  traveling  wave 
(A. 14).  When  x  =  x0,  the  transverse  part  of  the  electric  field  (3.36)  is 
the  electric  field  (3.38)  is  and  the  electric  field  (3.40)  is  When 

x  =  L?,  the  ratio  of  the  TE\o  voltage  associated  with  (3.38)  to  the  TE 10 
current  associated  with  (3.39)  is  Z2. 

Suitable  expressions  for  the  electric  and  magnetic  fields  due  to  the  M's 
of  (3.1)  and  (3.2)  are  given  in  the  previous  two  paragraphs.  We  are  nearly 
ready  to  substitute  these  M’s  and  their  magnetic  fields  into  (3.1)  and  (3.2). 
Letting  6  =  TM  in  (3.3)  and  using  (A. 4)  and  (A. 5),  we  obtain 

Ml™  =  A™(yV+)  (3.42) 

Letting  6  =  TE  in  (3.3)  and  using  (A. 16)  and  (A. 17),  we  obtain 

£E  =  a+,^)  (3.43) 
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Substitution  of  (3.42),  (3.43),  (3.31),  (3.33),  and  (3.35),  into  (3.1)  and  sub¬ 
sequent  application  of  the  orthogonality  (A. 26)  give 

Ipq 


Y  ' 


(3.44) 


Y 


1,ITE,ITM 


=  0 


y\,\TM,lTE  _  Q 
rij  —  u 


(3.45) 

(3.46) 


Y- 


=  < 


\,1TE,1TE 


i7io(c°sh  (7iq(Xi  -  x0))  +  ZtV;n^sinh(71o(Ii  -  x0)))6u 
w/i(sinh(7io(Ii  -  x0))  +  ZiY^  cosh  (tw(ZI  -  *«>))) 


.22 vahi 

CJ/i 


(p,?)  =  (1,0) 


,(p>9)  ^  (1,0) 


where  <5,j  is  the  Kronecker  delta  function: 


c  _  /  1,  *  =  J 
tj- \ 0, 


(3.47) 


(3.48) 


In  (3.44)  and  (3.47),  the  subscript  j  is  related  to  pq  as  in  (2.26).  The  subscript 
j  is  not  to  be  confused  with  the  other  j  in  (3.44)  and  (3.47).  This  other  j  is 

y/=l. 

In  a  development  similar  to  that  in  the  previous  paragraph,  we  obtain 


r2,2TM,2TM  _ 


jUtSy 


7pg 


r,2,2TE,2TM  _ 


Y: 


=  0 


y2,2TM,2TE  q 

*ij  U 


(3.49) 

(3.50) 

(3.51) 


y2,2TE,2TE 


J7io(cosh  (7io(^2  ~  Jo))  +  Z2^n~;  sinh  (710 {L2  -  x„)))6, 

w/i(sinh  (7io(L3  -  x0))  +  Z2Y{0E  cosh  (710(I2  -  *I)T) 

(?,?)  =  (  1,0) 

*  (ho) 


ill 


(3.52) 
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Chapter  4 

The  Admittance  Matrix  for 
the  Circular  Waveguide 

In  this  chapter,  Y?'a0nS  of  (2.30)  is  evaluated.  The  field  (£(3)(Q,M^), 
H^3\0.  Ml?))  is  radiated  by  the  magnetic  current  M^s  in  the  circular  wave¬ 
guide  with  the  apertures  A\  and  A2  closed  by  perfect  conductors,  with  the 
short  at  z  —  i3,  and  with  the  matched  load  at  the  other  end.  Otherwise 
stated,  is  the  field  that  would  exist  in  region  3 

of  Fig.  2.1  if  and  were  removed  and  if  MJ*  were  put 

where  —M^  was.  Our  first  objective  is  to  obtain  expressions  for  £(3)(0,  A£^) 
and  £(3,(Q,iO  An  expression  for  £*3)(Q,  A£p*)  is  needed  to  evaluate,  as 
indicated  in  the  last  two  sentences  of  Chapter  2,  the  field  in  the  circular 
waveguide. 

The  volume  density  of  magnetic  current  corresponding  to  the  surface 
density  M^q{<f>,z)  is  z)6(p  —  a)  where  6(p  —  a)  is  the  one-dimensional 

Dirac  delta  function.  This  volume  density  is  expressed  as 

M.p6q{<t>,z)S{p- a)  =  <i>,  z)  +  uzM^q(p,  <f>,  z)  (4.1) 

where 

«&(*>,  *>  =  (a*  '  Mifa,  *))«(*>  -  a)  (4.2) 

«£(/>,*.*)  =  -  a)  (4.3) 

Applying  the  Green’s  function  technique,  we  obtain 
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£(3,(0>MJ‘)  =  j\iz' ['  d#  [  R'iP'{M^,4,\z') 

•£(3)(W(!:  -  r'))  +  JO/,  2')£<3>(tt,  u,«(t  -  £))}  (4.4) 

££<3>(fl,AqJ)  =  J‘\  dz'  f  A#  £  z') 

•£'3,(Q,tt*<(E  -  r'))  +  r')£<3’(a,!i>(r  -  £))}  (4.5) 

A  (3) 

where  the  operator  R_  ’  differs  from  the  operator  R^  only  in  that  the  second 

« (3) 

argument  of  R  is  a  volume  density  instead  of  a  surface  density.  Similarly, 

«  (3) 

the  superscript  “A”  in  R  indicates  that  the  second  argument  is  a  volume 
density.  In  (4.4)  and  (4.5),  S(r  —  r')  is  the  three-dimensional  Dirac  delta 
function,  r  is  the  radius  vector  to  the  point  (p,  <f>,  z)  at  which  M^6 )  and 

H^lO.  M^sq)  are  evaluated,  and  r'  is  the  radius  vector  to  the  point  (p',  <p\  z'). 

The  fields  2£^(Q,  M^(r  —  L1))  and  j£^(Q,  t4(5(r  —  r '))  are  obtained  by 
adding  to  the  electromagnetic  field  of  (D.19)  and  (D.20)  the  reflection  due  to 
the  short  at  z  —  L *  The  electromagnetic  field  of  (D.19)  and  (D.20)  consists 
of  a  sum  of  fields  each  of  which  is  of  the  form  {R,  K)  where 

£  =  M*,V)£tan  +  n,£,)e-'l‘- ''  (4.6) 

£  =  (fltao  +  ».'(*.  (4.7) 

where  R and  2£tan  have  only  p  and  <f>  components.  Furthermore,  t (2,  z') 
is  given  by  (D.4),  and  7  is  either  7™  of  (B.24)  or  7jp£  of  (B.53).  In  this 
paragraph  and  in  the  next  three  paragraphs,  the  index  p  that  appears  in 
(D.19),  (D.20),  7™,  ard  7^£  is  not  to  be  confused  with  the  index  p  that 
appears  in  M^Sq,  A/Jp<?  and  M^q.  Similarly,  the  index  n  that  appears  in 
(D.19),  (D.20),  7£pM  and  7££  is  not  to  be  confused  with  the  index  n  that 
appears  in  The  reflected  field  due  to  (R,R)  is  a  field  (RT  ,H_r)  which  is 

proportional  to  (R,  R)z<z'.  Here,  the  subscript  z  <  z'  denotes  evaluation  at 
z  <  z' .  Adjusting  the  amplitude  of  ( Rr,H_T )  so  that  the  tangential  component 
of  R  +  El  vanishes  when  z  =  I3,  we  obtain 

Kr  =  (-R  un+uzEz)e^'-^] 

iT  =  (£tan  -  u,H2)e^+*'~2^ 
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(4-8) 

(4.9) 


The  sum  of  (4.6)  and  (4.8)  is 

£  +  £'  =  2c-'<i’-'>{£tan  sinh  (7(I,  -  *)) 

+UZEZ  cosh  (7(^3  -  2))},  z>z  (4-10) 

£  +  £'  =2e-’l1>-)coSh(7(i3-2')H-£tan+!4,£.}.  2  <  *'  (4.11) 

The  sum  of  (4.7)  and  (4.9)  is 

K  +  E  =  2e-<I>-'>{atancosh(7(£.3-r)) 

+uzHz  sinh  (7(1,3  -  z))},  z>  z  (4-12) 

R  +  Kr  =  2e-^-^cosh(7(l3-2,)){iitan-^^},  *  <  *'  (4-13) 

In  this  paragraph,  we  have  shown  that  the  short  at  z  =  L3  changes  the  field 
(22,  H)  that  would  exist  in  the  circular  waveguide  matched  at  both  ends  to 
the  field  (E  +  Er,K  +  Kr). 

Replacing  each  term  of  the  form  (4.6)  in  (D.19)  by  a  term  of  the  form 
(4.10)  or  (4.11),  we  obtain 

JE<3)(Q  ,i4*(i:-r')) 

1  -  ”  en(^)V;(fc^V)e-^M^-^ 

V  n=0p=l  XnpJn+  l(znp) 

■{-l±Mk™p)  cos  (n(0  -  <S>'))  sinh  (7™  (£3  -  *)) 

,  nMk™p)  sin  {n{4>  -  <£'))  sinh  {i™(L3  -  z)) 

+14  -jpT 

Knp  P 

,  k™  cos  (n{<t,  -  <t>'))cosh.(-r™(L3  -  z)) 


2  f-f.  (C)2fi^(Cl’')e^"(fa~,,)sinh(7r»£(^3  -  £» 
'  ^  ^  (*' 2  -  "2)(*LVR2(0 


n= 1 p=l 


n^n(fenp£p)cos(n(0-^))  T£ 


Kpep 


+  14 Jn(kn»p)  sin  (n(^  -  0'))}, 


z  >  z 


(4.14) 
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£(3|(W(n  -  e')) 

1  ^  £„ COsh  (7™(£3  -  O')) 

“  _  2-j  2-*,  _2  t2  /  \ 


n=0 p=l 


^-np^n+lC^np) 


{Mpi'(^npMp)  cos  (n(<£  -  <t>'))  -  % 


nMk™P)  sin  (n(<f>  -  <t>')) 


k™Mk™p)  cos  {n(4>  -  <t>')) 

+&* - — TM - } 


knp  P 


.  2  A  A  (k™)2nJn(k™p')e-<E^  cosh  (7^3  -  **)) 

*  £  «p2  -  n’XOUVnp) 


n^n(^p)cOs(n(<^-  <^'))  tI  /  i  TE 


kTJp 


~  ^J'n{klpP)  sin  (n(<£  -  <£'))}, 


2  <  2 


(4.15) 


Replacing  each  term  of  the  form  (4.7)  in  (D.20)  by  a  term  of  the  form  (4.12) 
or  (4.13),  we  obtain 

£(3,(W(E-t')) 

_  y.  c.(g)2  cosh  (7™(£.  -  --)) 

^n=Op=l 

,  ..  ./ / iTW . _ /_/  . 


+  14-7n(inpWp)cOs(n(0  -<?'))} 


J2  ^  ^ 

^  ntl  Ptt  (x;p2  -  n2)(klEp')J2n(x'np) 

f/  ti  fiTE  ^  ,  nJn(^£p)cos(n((?  -  0')i 

■{(upvfn(fcnp  />)  sin  (n(<£  -  (i> ))  +  14 - 2 — jjg - - - 

knp  P 

•  cosh  (7JpE(I3  -  2 )) 

klpMknpP )  sin  W  ~  <£'))  sinh  (tJp£(^3  -  z)) , 


},*>z'  (4.16) 
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fl<3)(W(n  -  tf) 

-jk\ g  g  *„ cosh(-,™(i3  -  z')) 


7' 


TA/  2  r2 
np  ^np^n+l 


(*np) 


r  nJn(fcnp  p)sin(n(0-^))  r,,f™  v  ,  ,, 

•K - F  lrTM  , -  +  n*J»(fcnp  />)  COS  («(<£-<£))} 

Tnp 


ki‘“p 


’np 
s, T£/ 


,  j2  g-  g  7y(QW.(tj[.V)e-g(t»-»)  cosh  (7™(L3  -  z')) 


"P 
J  2 


(*np2  -  "2)(OVn(0 


1fUJPn=  1  p=l  \~np  "  /Vnp  r  /"nV-npy 

•{Mp4(^np£/>)  sin  («(<£  “  <£')) 


nJn(knvP)  cos  (n(0  -  0')) 


+14- 


kTE. 


’np  P 

knpMknpP)  ^ n(w(^  -  (ft'))  ^  , 

2  ^yTE  J  ’  2  ^  ^ 

/np 


(4.17) 


The  fields  £  (Q,Ji*6(i:-  t'))  and  H.  (Q,li*6(E.  ~  r'))  are  obtained  by 
adding  to  the  electromagnetic  field  of  (D.55)  and  (D.58)  the  reflection  due 
to  the  short  at  z  —  £3.  Aside  from  the  8{z  —  z ')  term  in  (D.58),  the  electro¬ 
magnetic  field  of  (D.55)  and  (D.58)  consists  of  a  sum  of  fields,  each  of  which 
is  of  the  form  (RTE ,  ifT£:)  where 


KTE  = 

KTE  =  («(*,*')£&  +  uzH™)e-W -*'1 


(4.18) 

(4.19) 


rn  ri  p 

where  £tan,  and  iftan  bave  only  p  and  4>  components.  The  reflected  field 
due  to  {E^ ,HJE)  is  called  (£T£r,  KTEr)  and  is  given  by 


ETEr  =  -£^e7"',E(z+z'-2£'3) 

RTEr  =  (Hj£-  u,HjE)e-’"i‘*1'-2L>) 

The  sum  of  (4.18)  and  (4.20)  is 


(4.20) 

(4.21) 


ETB+ETEr  =  2£S,c-*E(‘>-')sinh(7jpE(r3  -  z)),z  >  (4.22) 

Ere  +  HTE’  =  2£gile-'”<t»->  sinh  hi®(£ 3  -  z')),  z  <  z'  (4.23) 
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(4.24) 


The  sum  of  (4.19)  and  (4.21)  is 

KTE  +  gTEr  _  2e-,.V(ta-.') cosh  {7TK(£3  _  2)) 

+U,HJB  sink  ('tlfiU  -  z))},  2  >  2' 

f£TE  +  i£TEr  =  2e-’W‘->  sinh  (7” (I,  -  =')) 

{-flS,  +  2  <2'  (4.25) 

The  short  at  2  =  jL3  changes  the  field  {RTE  ,KTE)  that  would  exist  in  the 
circular  waveguide  matched  at  both  ends  to  the  field  (£.TE  +  RTEt  ,Kte  + 

KTEr). 

Replacing  each  term  of  the  form  (4.18)  in  (D.55)  by  a  term  of  the  form 
(4.22)  or  (4.23),  we  obtain 

£<3|(Q,u,%-r')) 

=  1  £.  g  e. (^f£,(kTJ sinh  (7 -  2)) 


n=0 p=l 


nJn(knE  p)  sin  (n(<f>  4>))  t/hTE  \  >  ti  nwt 

U±p - ZteZ - +  *4 >JniknpP )  COS  -  <t> ))}, 


kTEi 
Knp  / 


z  >  z 


(4.26) 


1,(3) 


'(Wfe-cf)) 

=  1  g  g  fn(C ££n(*iV)e-,”(I,-*)  Sinh  (7Jf  (*3  -  *')) 


n=0  p=s  1 


JTE(  >  2  „25  t2/  '  \ 
7np  (^np  n  )c'n(a'np) 


nJn{k  p)  sin  (n(0  —  <j> ))  ti/iTE  \  r  ti  i>\  m 
•{mp - 2 — - +  ^t,Jn{knp  p)  cos  (n(<£  -  <t> ))}, 


*?/» 


<  z 


(4.27) 


Replacing  each  term  of  the  form  (4.19)  in  (D.58)  by  a  term  of  the  form  (4.24) 
or  (4.25)  and  retaining  the  S(z  —  z ')  term  in  (D.58),  we  obtain 


tn{klv?  Jn{klE  p‘) 


fr(3)/n  r/  in  J  O  V"  Cn  V  Knp  / 

k  da.  (c  c.))  -  TUfl  E  E  (<i  _  j 


nz8(z  -  z\ 


,Mknpp)  cos {n{<t>  -  <t>')) 


klE 


+ 


H+,z'  >  z 
H_~ ,  z'  <  2 


(4.28) 
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where 


u-  If  T/ / 1 te  \  ,  ,i  i/\\  nJn{k™p)sin(n{4>-4>')) 

a  -<  UlMKvP)  cos  -<£))-  JU -  kTEp - 

•  (£3  -  z)) 

kl pMknpP)  cos(n{<t>  -  #))  sinh(7jp£(  C3  -  z))  1 


(4.29) 


it+  =  [v.ttgpl-W*  -  *0)  (*  ~  ^ 

_^g^g£)co!W-^|  sinh(7rE(i3  _  0)  ,,.30) 

7np  J 

Before  substituting  (4.14)-(4.17)  and  (4.26)-(4.28)  into  (4.4)  and  (4.5), 
we  obtain  expressions  for  and  M?£q.  Substitution  of  (A. 10)  into  (2.13) 
and  (2.14)  gives 


1  nTM  / 1  ^  2 IT  {  q  .  ,p*y'1+ s  /97rz+N 

^  (*'z)  =  -^ia<csm(— )C0S(  — ) 

/  ,v* ,P\,suLt0'  ,pn ry'T+,  .  ,qirz+^\ 

-Uz(-l)  (j)("T“)  cos(— ^— )  sin(  -  ) 


(4.31) 


where  7  is  either  1  or  2.  Moreover,  y7+  and  z+  axe  given  by  (2.15)— (2.17). 
Finally,  from  (A. 8),  kpq  is  given  by 


W<  f)J  +  <T)! 


(4.32) 


Substitution  of  (A. 23)  into  (2.13)  and  (2.14)  yields 

?sin(^)cos(^ 


&p,  V  46c  (  6 


,  /  tvr/?wsin^»N  ^  moo, 

+Uz(-1)  (-)(—: —  )cos(— — -)sin( - )  >  (4.33) 

c  <Po  0  c  J 

where  ep  is  Neumann’s  number  given  by  (A. 21).  Substituting  (4.31)  into 
(4.2)  and  (4.3),  we  obtain 


^™(p,  <M)  =  ^=(-)sin(^|—  )cos(^-)<5(p  -  a) 

kva  vbc  c  b  c 


(4.34) 
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i *~<tm t  i  \  2?r(-lk  p  sin0o,  /p7ry7+N  .  ,q^  xc/  ,  ,A  OCN 

4>,z)  =  — y=-(£)(-^—  )cos(— —  )sm(— —  )8{p-a)  (4.35) 

Substituting  (4.33)  into  (4.2)  and  (4.3),  we  obtain 

A I£,EM,Z)  =  -■pyS(£)sin(mZl)cos(2^)^-a)  (4.36) 


pny 


n+ 


•sin(- - )6(p  —  a)  (4.37) 

c 

Substituting  (4.34)  and  (4.35)  into  (4.4)  and  performing  the  integration 
with  respect  to  p\  we  obtain 


{ 


/+  r 


/9n  .  ,9**  X 

( - )  sm(  -  )  c°s(  -  ) 

CO  c 


£|3|(Q,“*A(t  -  r')) 


,  ,\-,fP\/sin <f>°\  ,V*y  s  ■  ,95rz'+J^(3)/n  /u 

-(-1)  )cos( — —  )sm{— — )  |£  (Q,uz<5(r-r  )) 

where  7  is  either  1  or  2  and,  recalling  (2. 15)— (2. 1 7) , 

y'1+  =  (7T  -  <j>')x0  +  f 

y'2+  =  ^xo  +  l 
/+  - 


p  —a , 


z-  =  z'  + 


01  =  (2  -  7)?r  -  0O 

02  =  (2  -  7k  +  0o 

Similar  substitution  of  (4.36)  and  (4.37)  into  (4.4)  gives 


(4.38) 


(4.39) 

(4.40) 

(4.41) 

(4.42) 

(4.43) 


£,3,«E> = 


(-)sm( — - — )cos(— ) 


,(3) 


K  '((WU-n')) 


,'i+ 


,  /  ,w,9wsin<^  ,p*y"' ,  .  ,9*^'  , 

+(-l)  (-)(— 7—  )c°s( — 7 — )sin( ) 

C  <po  0  C 


KW(Q,llAr-r!)) 


(4.44) 


J  p'  =  a 
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Since  (4.5)  is  (4.4)  with  R^  and  E replaced  by  and  respectively, 

the  equations  obtained  by  replacing  R ^  and  R ^  by  and  2£(3)  in  (4.38) 
and  (4.44)  are  valid: 


{ 


i£(3)<a ,<■") = --^=  iz  r  i* 

kpqy/bc  J-\  J<tn 


•  ,p*y,y+\  ,<l‘K2'Jr^ 

(-)sm( — - — )cos( - ) 

CO  c 


K<3>U1,u<,6(l  -  r'))l 


J  o'  — 


p'-a 


-(-1)7(?)(-T^£)cos(^-T — )sin(^— )  K^i&iiJir-r'))  \  (4.45) 

0  C  L  p'—a  | 


4>o 


{ 


r'+ 


,Px  •  fP*y  X  ,q*z-  x 
(-)sm( — - — )cos(  — — ) 


a(3,(a,M(r-r')) 


p'=a 


+(-l)7(-)(~T—  )cos( — 7 — )  sm( - ) 

c  <t>0  o  C 


h!3\h,u,6(l  -  c')) 


(4.46) 


Replacing  (n,p)  in  (4.14),  (4.15),  (4.26),  and  (4.27)  by  (r,  s)  and  then 
substituting  the  resulting  expressions  into  (4.38),  we  obtain 


nTE<t> 


kpqVbc\c  rt£  fr{  2 xj,J?+l{xr3) 
{kTEyrjr{kTEa)K 

c  fr[  ftj  {kJE a)(x'r1  -  r2)J}{x'T,) 

p  sin^,  ”  ”  e1(TO(e^ ) 

6  to  27j,£«,2  -r2)Jr2(x/j  J 


where 


fi  ,g7rz'+ 

=  Lcoa(-r> 


ppTMib 

R?te*  =  J*  cos( 

£7T£*  =  sjn( 


2 


g7TZ 


'+ 


(?7rz 


'+ 


£,tm*+,z,  >  2 
<  2 

gyTE<t>+,  z,  >  z 

r:te*-,z'  <  2 
w  R?TEz+,z'>z  ' 
\  EPTEz-,zi  <  z 


-) 


Jdz' 

>  dz' 
dz' 


(4.47) 

(4.48) 

(4.49) 

(4.50) 
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In  (4.48)— (4.50), 


£PTm~  =  -e^.M^3-z‘)  |(j LpF2J'r{k™p)-u 


<PxrJr{k™p), 


■s\nh(i™(L3  -  z))  -u. 


wr\  rs  rj  *4  fcTM p 

^kJMjrjkJMp)  cosh(7™(L3  -  z)) 


IT™**  =  cosh (7™(L3  - 

^'rjr(i™p)  .  p'k™Jr(kZ»p)\ 

~**~k JFT~  +  U’ - IfF - 1 

grTE*-  _  _e-T&(L,-»’)  sinh(7rT/(Z3  -  2)) 

E?te++  =  cosh(7rT/(£3  -  2'))e"%T*£(i3“2) 

/  ^rJrjkJfp)  x  ,  TE  1 

kTEp  ^ri^ra  p)  J 

^T£,-^c^J4(n,-I')3i|ih(7TS(Za_^(_i)7  .  _ 

£PTEz+  =  sinh(7 ™{Lz  -  z'))e-^E(Li-z\-\y 


where 


+  ^J'r(klEp) 


<fx  =  f  sin{^~ — )  sin(r(<£  —  <f>'))  dtf 
J<bi  b 

4?2=  [**  sin(^—~—)cos(r(4>  —  <p'))  d<t>' 

U\  b 

4>l3  =  /  cos(  — — 7 — )sin (r(<j>  —  <t>'))  d<t>' 
J<b  i  0 

4P*  =  /  cos(^~ — )  cos(r(<t>  —  &))  d<f>' 

J<t>  i  o 


(4.51) 


(4.52) 


(4.53) 


(4.54) 


(4.55) 


(4.56) 


(4.58) 


(4.59) 


(4.60) 


Replacing  (n,p)  in  (4.14),  (4.15),  (4.26),  and  (4.27)  by  (r,  5)  and  then  sub- 
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stituting  the  resulting  expressions  into  (4.44),  we  obtain 

nf3l/n  tryTE\  _  4a  [Zfa  (  P  ^  ^  tr(kJsM)2  _ * 


*-w>~gi/S{tgs 


+?EE 


2x2jJr2+1(xr#) 

(k™)*rJr(k™a)E:TE* 


b  ^AV.E«)(Xrl  ~  r2)J?(x'ra) 
q  sin^  ”  ”  er{kJEfJr{k^g)ErT^\ 
+lcn  <*0  ZlTE(xrl  ~  r2)y2(x'j)  / 


(4.61) 


Replacement  of  (n,p)  in  (4.16),  (4.17),  and  (4.28)  by  (r,s)  and  subsequent 
substitution  of  the  resulting  expressions  into  (4.45)  lead  to 


fj( 3W0  fupTM\  _  4Ja  J  V'  V  er(kJaM)2J'r{kJaMg)irT 

K  ~  UJflkp^y/bc  \  C  ££  2 Wz’^ixr.) 

q  ^  ^  7^  (£*)Wr(£MfTrE* 

c  £i  (kZEa)(x'r2  -  r2)J2{x'ra) 

i  (  ir(pHsin^i,r y* er(kjaE)3jT{kj3Ea)irTEz \  ,4  62x 


rsO  1=1 


where 


ITte*  =  sin(^— — — )  [iT7 

7-  i  c 


/  JF™++,z!  >  z 
c  )\TTm~^<z 
qirz,+ .  /  jrTE*+,Z'>z 


■yTEt± 


irTEz+,z'  >  z 
irTEz-,Z'  <  z 


(4.63) 


(4.64) 


(4.65) 


In  (4.63)-(4.65), 

2T™*'  _  e-'r”M(L3-z')cosh(7rr;W(X3 -2)) 

=  cosh(7™(L3  ~  z'))e-^M(Li~z) 

\  <jPlrJr{k™p)  _  r'n.TM.x 

•|^p - - +  “♦*  Jr(^,  />) 


(4.66) 


(4.67) 


H.™-  =  | + 1 

,  ,  te,t  xx  Jr{kJsE p)  sinh(7^£’(X3  —  z)) 

•  cosh(7;4  {L3  -  z))  -  Ux - “tk - 

ITS 

££*te*+  =  cosh(7rT/(i3  -  z'))t-^E(U-z] 

L  df1  J'  ( kTE  o)  +  u  .  u  flSJASA 

'  | U-PV  Jr\Ka  Pj+24  kTEp  JTE 


(4.68) 


(4.69) 


if 'TEz±  =  «,*(*  -  ^)— ^14^"  *  (4-70) 

irT£-  = 

,  ,  tb/  r  xx  (k?sEp)  sinh(7rr/(L3  -  2))  ] 

■  cosh(7„  (Z-3  -  2))  -  u* - ^ - f  (4.71) 

/rs  J 

if  rE,+  =  -sinh(7rr/(l3  -  2'))e"7?/(L3~*) 

•  - 1 ■■)  (4.72) 

Replacing  (n,p)  in  (4.16),  (4.17),  and  (4.28)  by  (r,s)  and  then  substituting 
the  resulting  expressions  into  (4.46),  we  obtain 

lf(3)(r)  xf-yTEx  __  4ja  ltptq  j  k2p  ^  y,  e T(k™ )2 J'r(k™ a)IP 

P  v~-  f-  7 lE(k™YrJ,(k™a)S:TE* 
T>kk  (k™a)(x'*-r*)J}(x'r,) 

-(-‘)  (-)(—) EE  2(l,?_rW,(i;<)  }  (4-'4> 

Substitution  of  (4.51)-(4.56)  into  (4.48)-(4.50)  gives 
e™*  =  {-r™>  Sinh(7™(L3  -  z))  + 

+  {2™>  cosh(7™(i3  -  *))  +  ^-*klMUk™P)  (4  74) 
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£pTE4>  _  {_ZTE1  sinh('y^£'(Z,3  -  2 ))  +  zT  E2  e~'1',E 

f  ^rJrjkJfp)  _  !  J,,,TE-S 

'  kTE p  ^r\Ks  P) 

g,TEz  _  {^3  sinh(7T£^3  _  .,))  +  zTE4e^B{L3. 

(  lY*fu  |  y-»4  V(i-TE  \ 

•(-!)  jli, - "kj^~p - +  ^  J^kra  P 


where,  for  6  equal  to  TM  or  TE, 


zsx  =0,  2  <  -i 


zSl  =  f  e"^(Li'z,)  cos(~~~~) <fz,)  -f  <  2  <  f 
y-|  c 

251  =  ^  cos(^—  )dz',  z  >  2 

z52  =  1 2  co3h(7^(I3  -  2'))  cos(^i_)  dz\  z  <-  § 
y-f  c 

»°  =  /*  cosh -  *0)  c°s(— )  <fe',  -f  <  *  <  f 

«/X  C 


2"  =  0,  2  >  f 


2TE3  _  2  <  -  § 


2  <  — ; 


zrE3  =  j\  e-3”'1’-'*  sin(  ^ )  dz',  -  f  <  z  <  f 

zTE 3  =  J_\  e->.W.-.')  sin(^)  rfz',  2  >  | 

zTE' =  J’ ^sinh(-rj£(£3-z'))sin(^—  )dz\  z  <  -f 

zTEi  =  l2  smh{iJE(Lz  -  z'))sin(^—)iz',  -|  <  z  <  | 
y*  c 

,r£4  _  n  _  ^  c 

2  =  U,  2  >  j 

Substitution  of  (4.66)-(4.72)  into  (4.63)-(4.65)  gives 
££iTM<t>  =  |2TM1  cosh(7TM(L3  _  z))  +  2m2e-^(L3-„J 


(4.75) 


(4.76) 


(4.77) 

(4.78) 


(4.79) 


(4.80) 

(4.81) 

(4.82) 


(4.83) 

(4.84) 


(4.S5) 


(4.86) 


(4.87) 

(4.SS) 
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fT™*  _  |2T£1  cogh(7TE(l3  _  +  2TKe-%r.£(L3-*)j 


(4.89) 


L.e'j'' 


’TE  \  ,  ..  ^2rJT(kJ3Ep) 


(KsEp)  +  u*- 


klEp 


{-tTa  srah(-,J£(I3  -  *))  +  ^«e-,5*(L,-,)} (4.90) 

IT™*  =  {z™3  COsh(~flE(L3  -  z))  -  zTE4e-^?,B(L3-z)} 

•  { M’mTv)  -  + {*T£5 

-3T“sinh(7™(I3  -  *))  -  z™e-*’«>-‘>}ufik™JTf™p)  (4.91) 


where 


_TES  _ 

Z 


T?  <&*{*?■) 

0, 


,  -f<4< 
otherwise 


(4.92) 


In  (4.92),  z+  is  given  by  (2.17). 

So  far  in  Chapter  4,  we  have  found  that  £l(3)(0.  £(3)(0. 

and  H}3)(b,M£E)  are  given  by  (4.47),  (4.61),  (4.62),  and 
(4.73),  respectively.  In  (4.47)  and  (4.61),  Jj?™*,  £n'r£'<*>,  an(j  j^tez  are 
given  by  (4.74),  (4.75),  and  (4.76),  respectively.  In  (4.62)  and  (4.73),  fp™*, 
H?TE4>,  and  JpTEz  are  given  by  (4.89),  (4.90),  and  (4.91),  respectively. 

The  quantities  <jPl,  <fp2,  <p3,  and  <pA  which  appear  in  (4.74)-(4.76)  and 
(4.89)-(4.91)  are  evaluated  in  Appendix  E.  These  quantities  are  given  by 
(E.10)-(E.13)  in  which  <j>^\  <j>^\  <^3\  and  are  given  by  (E.23)-(E.26). 

The  quantities  z™1,  zT™2,  zzEl,  zTE 2,  zTE 3,  and  zTEA  which  appear 
in  (4.74)-(4.76)  and  (4.89)-(4.91)  are  evaluated  for  (  — |  <  z  <  |)  in  Ap¬ 
pendix  F.  As  indicated  in  Table  F.l,  the  quantities  z™1  and  z™ 2  are  suit¬ 
ably  given  by  (F.25)  and  (F.26)  when  7™  is  purely  imaginary  and  by  (F.32) 
and  (F.33)  when  7™  is  purely  real.  In  (F.25)  and  (F.26),  0&T,  is  given  by 
(F.27),  q5~  is  given  by  (F.23),  and  qs+  by  (F.24).  The  quantities  zTEl,  zTE2. 
zTE3,  and  zTEA  are  suitably  given  by  (F.25),  (F.26),  (F.34),  and  (F.35)  when 
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lJE  is  purely  imaginary  and  by  (P.32),  (F.33),  (F.41),  and  (F.42)  when  ~fJE 
is  purely  real. 

If  z  <  —  then  z™1  =  zTEl  =  zTE 3  =  0  while  z™2,  zTE 2,  and  zT£4  are 
given  by- their  expressions  in  Appendix  F  with  z+  replaced  by  0.  If  z  >  |, 
then  z™2  =  z7^2  =  zTE4  =  0  while  z™1,  zTEl ,  and  z7"^3  are  given  by  their 
expressions  in  Appendix  F  with  z+  replaced  by  c. 

In  view  of  (4.89)-(4.91),  substitution  of  (4.31)  and  (4.62)  into  (2.30)  gives 


y3,aTM,-lTM  =  T  |2|(5,  -  5,)  +  ~S3 


(4.93) 


where 

T 

51 

52 

53 

54 

55 


8irja2 

kjnnk^Ufibc 

-  «  er(k™)2J'r2(k™a)z^2 

2i™x?,j?+  i(*«) 

^  7 ™(kZEyr2J?(k™a)z(2)r-'2 

hh  WEa)2{x'2-T2)J2{x'ra) 

sin^  f,  (■ kZE)3rJ?(k™a)z 

K  )K  0»  (fc"a)(li2_r»)J?(x;.) 

sin^0  cr(^)3rJ2(^a)zt4)^3 

1  j  1  2(k™a)(x'2-r*)J2(x>s) 

,  na^/ sin  ^0X2  f  f  er{kJsE)4J2{kEEa)z{5)<j>a'1‘i 
1  j  1  00  ;  ££s  27«(*,„a-r»)J?(*;,) 


(4.94) 

(4.95) 

(4.96) 
(4-97) 

(4.98) 

(4.99) 


In  (4.95)-(4.99), 


r*ml 

,,i  /m7ry“+ 

<fp  cos( — - — )d<p 

^3  0 

(4.100) 

=  i 

j-y2  •  /"*XSf<H\ 
p7  sin( — - — )  d<p 

^3  b 

(4.101) 

(m7“+)^ 

^3  0 

(4.102) 

^a^4  =  y 

tp*  cos( - ; - )  dtp 

i3  0 

(4.103) 
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where 


(4.104) 

(4.105) 


<t>3  =  (2  -  Q)5T  -  <j>0 
<f>  4  =  (2  —  a)ir  +  <f> o 

Still  in  (4.95-(4.99), 

zW  =  P  { z Tm  cosh(7™(I3  -  *))  + 

j-t  < 

mrz+ 

•cos( - )  dz  (4.106) 

c 

*'2'  =  J\  K*1  -  *))  +  zTE2e-’'r-‘<l>-)} 

^  3 

•  cos( U— — j  <fz  (4.107) 

c 

z^  =  /2  {— zTE1  sinh(7^£'(X3  —  z))  +  zT£2e-'|,"^£'3-** j 

■sin(-— ---  )dz  (4.108) 

c 

z<4>  =  j\  {*r,Oc08h(7™(i,  -  *))  -  *«•«-£*(*»-*>} 

J  7 

717TZ  + 

•cos( - )dz  (4.109) 

c 

z<!'  =  /S  {;rES  -  zr“«ah(7;B(Is  -  *))  -  *r*4e-,£«(£,— )} 

-2  1 

•  sin( - )  dz  (4.110) 

c 


In  view  of  (4.89)-(4.91),  substitution  of  (4.33)  and  (4.62)  into  (2.30)  gives 


y3,aT E  ,*yT  M 


(Si  -  Si)  -  ^S,  - 

c* 


mP  c  ,  nP  n 
b*  64  +  6c 


(4.111) 


In  view  of  (4.89)— (4.91 ),  substitution  of  (4.31)  and  (4.73)  into  (2.30)  gives 


■y3,aTM,-rTE 


2£,c  c  ^  ,  ailq  ,  Ulq  ■ 
6c(5i-52)+  62  53  +  c254  + 


Similar  substitution  of  (4.33)  and  (4.73)  into  (2.30)  gives 
y^TE.TB  =  _  S,)  -  ^ S 3  +  ^ S ,  -  ^S,}  (4.113) 

In  the  previous  paragraph,  we  found  Y?'a™'y™ ,  y3,aTE,-,TM , 
y3,aTM,-fTE ,  ^  y3,aTE,yTE  tQ  bg  gj^  by  (4.93),  (4.H1),  (4.112),  and 

(4.113),  respectively.  In  these  equations,  T,  Si,  52,  S3,  S4,  and  S$  are  given 
by  (4.95)-(4.99).  The  quantities  <j>a'lX ,  <f>ay2,  <t>ayZ,  and  <^an'4  that  appear  in 
(4.95)-(4.99)  are  evaluated  in  Appendix  E.  These  quantities  are  given  by 
(E.3l)-(E.34)  in  which  <fy\  <f>^2\  <j>^\  and  are  given  by  (E.23)-(E.26). 
Also  in  (E.31)-(E.34),  <£al'Yl,  4>a2yl,  <t>aly2,  and  <f>a2y2  are  given  by  (E.46)- 
(E.49)  and  (E.53)-(E.56).  The  quantities  z^\  z*2),  z^\  and  z^  that 

appear  in  (4.95)-(4.99)  are  evaluated  in  Appendix  F.  These  quantities  are 
given  by  equations  whose  numbers  are  listed  in  Table  F.l. 
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Chapter  5 

The  Excitation  Vector 


In  this  chapter,  1°^  of  (2.24)  is  evaluated.  Approximate  expressions  for 
i£(3)(,Z.imP,Q)  and  have  to  be  found  before  the  integral  in  (2.24)  can  be 
evaluated. 

In  (2.24),  is  the  magnetic  field  due  to  radiating  in  the 

circular  waveguide  with  the  apertures  A\  and  A2  closed  by  perfect  conductors, 
with  a  perfectly  conducting  wall  at  z  =  and  with  a  matched  load  at 
the  other  end.  If  the  apertures  were  present,  the  only  2  traveling  wave  at 
2  =  —  |  contained  in  the  field  of  iimP  would  be  the  unit  amplitude  2  traveling 
TMqi  wave.  Since  closing  the  apertures  produces  no  2  traveling  waves  in  the 
region  for  which  2  <  — the  only  2  traveling  wave  contained  in  the  field 

at  z  —  —  |  is  the  unit  amplitude  2  traveling 
TMqi  wave.  The  field  of  this  wave  is  ( ,  H%™e+ )  given  by  (B.l): 


pT  Me+  _ 

£-01  ” 

rrTMe 4-  _ 

2*01  ~ 


Z0™eoe™e(p,  +  u; 


(fc™)20™e(l,0)e  100 

jut 


tm7 
1  z 


(5.1) 

(5.2) 


Here,  k^M  and  ip™e{p,  </>)  are  given  by  (B.7),  eJiWe(/5,  <f>)  by  (B.22),  h ™e{p.  6) 
by  (B.23),  and  by  (B.25)  in  which  7™  is  where 


0™  =  y/k*  -  (k™?  (5,3) 

Now,  f)™  is  purely  real  because  it  was  assumed  that  the  T.V/0I  mode  prop¬ 
agates  in  the  circular  waveguide. 
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The  field  (£(3)(^imp, Q), it(3)(iZimp, Q))  is  (i£,Me+ , i££Me+ )  of  (5.1)  and 
(5.2)  plus  the  reflection  of  (£j1We+,i£^1Me+)  from  the  conducting  wall  at 
Z-3.  Taking  the  reflected  field  proportional  to  (E™f~ .  H™e~ )  of  (B.2)  and 
requiring  the  p  component  of  JE^(iZimp,Q)  to  vanish  at  z  =  L$,  we  obtain 


E?'U toP.fl)  =  {7™sin(/30™(i3  -  *))&"*(/>,*) 

u;c  *• 


4>)co s($f(L3  -  z))} 

(5.4) 

i£(3)Uimp,fl) 

=  cos (S™ (Lz  -  z))lSf'M) 

(5.5) 

Equation  (B.25)  was 

used  iD  obtaining  (5.4).  From  (B.7),  we  have 

j^TAfe/  i\  _  MW. P) 

(5.6) 

where 

lTM  _  *01 
kQ1  -  a 

(5.7) 

Applying  [4,  eq.  (D-15)]  to  (B.22)  and  (B.23),  we  obtain 

. TMr /  i\  MW  P) 

(5.8) 

fjTMe,  J\{k™  P) 

(5.9) 

Substitution  of  (5.6),  (5.8),  and  (5.9)  into  (5.4)  and  (5.5)  gives 


O  -j0™L3 

£<3'Uin,P,Q)  =  {a»7o,;"^(io™^)=in(/J0™(i3  -  .-)) 

Jo(k™p) cos (/3™(I3  -  =))}  (5.10) 

a|3|uimp,!i)  =  ^;X(~x^)!4~'1(t°™  ^  c°3(<?™(i'3  ~ z])  (5-u) 

Although  not  needed  to  evaluate  expression  (5.10)  is  needed  to  evaluate, 
as  indicated  in  the  last  two  sentences  of  Chapter  2,  the  field  in  the  circular 
waveguide. 


_ j 
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Replacing  (7 ,p,q)  by  (a,  m,n)  in  (4.31)-(4.33),  we  have 


AC™(<M)  = 


2  7T 


kmn'/bc 


{n  .  .m7rya+.  .  uttz* 


,»+ 


./  iwmwsin<^x  ,mTy'*~r  .  ,nnz-  sl  inN 
— 3A«(— 1)  )cos( — - — )sm(— —  )>  (o.l2) 


2?T  /  f  r7T 


.  rrnry  +  ,nirz+ 

sin( - - — )  cos( - ) 

0  c 


,a+ 


/  ,  x«/nv/Sin  c6ON  .  m7ry“",\  .  n7rzn  ,  ,  / 

+iiI(-l)0'(-)(-^)Cos(— )sin(  — )|  (5.13) 


where 


(5.14) 


In  (5.12)  and  (5.13),  ya+  and  z+  are  given  by  (2.15)-(2.17). 

Substituting  (5.11)— (5.13)  into  (2.24)  and  integrating  with  respect  to  yQ+ 
and  z+  rather  than  <f>  and  2,  we  axrive  at 


aTM  8<ft0n  a2TT 


ralM  _ 

ii 


-i0lML3 


kmn  be 


VsmZcaxe  1  01 


jaTE  Wam  c  /a27T£me 
’  *mn«VV  46c 


ny,mZcc ne  ^ 


(5.15) 

(5.16) 


where 


In  (5.18), 


1 

J/sro  —  ^ 

r*=a=>*~ 

(5.17) 

_  1 
^ccn  — 

C  . 

r  cos(f3™( 

JO 

L3  —  2+))cos( - )  dz  + 

c 

(5.18) 

Lt  = 

r  c 

=  Z,3+2 

(5.19) 

(5.15)  and 

(5.16),  we  used  (5.7)  and 

Xo 

6 

"  20o 

(5.20) 
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If  mt  ±  0™c  =  0,  then  the  right-hand  side  of  (5.23)  must  be  replaced  by  its 
limit  as  nir  ±  approaches  zero.  To  render  the  value  of  this  limit  obvious 
and  to  avoid  roundoff  error  when  j nir  ±  0™c\  is  small,  we  recast  (5.23)  as 


_  sin(n7r  -  /%i Mc)  cos {P™  L$ )  -  2  sin2( sin(/3j1Af L3  ) 

2(n7r  -  0™c) 

.  sin(nff +  ^1Mc)cos(/3^MIJ)-l-2sin2(2^fu-^)sin(^1MIJ)  ,B 

+  2(n7T  +  0™c)  (  j 


Zccn  — 


Now,  If™  and  IfTE  are  given  by  (5.15)  and  (5.16)  with  ysm  amd  zccn  given 
by  (5.21)  and  (5.24). 
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Appendix  A 

Modes  of  the  Rectangular 
Waveguide 


Consider  the  rectangular  waveguide  whose  cross  section  is  shown  in  Fig.  A.l. 
In  Fig.  A.l, 


V+  *  V  +  |  1 

*+«*  +  §  J 


(A.l) 


where  y  and  z  are  shown  in  Fig.  1.1.  Four  kinds  of  mode  fields  can  ex¬ 
ist  in  the  waveguide  of  Fig.  A.l.  These  mode  fields  are  (£^f+,  2£^f+), 
(j Elf+,2£^+),  and  (£ZSr,l£Sr).  Here,  £  is  the  electric 
field  and  E_  is  the  magnetic  field.  The  superscript  “T M"  denotes  transverse 
magnetic,  UT En  denotes  transverse  electric,  “-fi”  indicates  that  the  wave 
travels  in  the  -fix  direction,  and  ”  indicates  that  the  wave  travels  in  the 
—x  direction.  Here,  x  is  the  rectangular  coordinate  measured  in  the  u.y  x 
direction  where  u.y  and  uz  are  the  unit  vectors  in  the  y  and  z  directions, 
respectively. 

From  the  analysis  in  (4,  sec.  8-1],  we  obtain 


EZH*  =  z™£X(y+,z*)e-""x  +  ur 

Sl"+  =  /^n(y*,z+)e-^- 


(A. 2) 
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Figure  A.l:  Cross  section  of  the  rectangular  waveguide. 


and 


vV™  +  =  o 


TM 


subject  to  the  boundary  condition 


1.2  T+\p~lmnX 

=  -z™£X{y\*+Y"'nX  +  uMmnU---)e 

jut 

where  Mx  13  the  unit  vector  in  the  x  direction.  Moreover, 

£^(!/+,z+)  =  -Y0™(j,+,Z+) 

k™(v\z+)  =  -u,  X  Y'/>,™(!/+,--+) 

The  wave  function  V>™(y+,  2+)  satisfies 


(A-3) 


(A. 4) 
(A.r>) 


(A. 6) 


V>™  =0 


(A. 7) 
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on  the  walls  of  the  waveguide.  Solutions  to  (A. 6)  and  (A. 7)  are 


kin  =  (™)2  +  (™)2 

0™(y+,*+)  =  T~  2  /7-  sip(  )  sin(  ) 

KmnVbC  t>  C 


m  =  1,2,  •  •  • 
n  =  1,2,-- 


(A. 8) 


The  preceding  ip™  's  normalized  so  that 


j£V  />+  {E<£5V,*+)}  ■  {SiOs/V4-)}  =  1  (A.9) 

Substituting  °f  (A. 8)  into  (A. 4)  and  (A. 5)  and  taking  £  with  respect 
to  the  cordinates  y+  and  z+,  we  obtain 


^(yV+)  = 


2ir 


kmn  Vbc 


{“■ 


m 


cos( 


m7ry^ 


')  sin( 


nirz1 


+li,-sin(- 

C  0 


-)  cos(- 


A^(i/V+)  = 


2tt 


{n  .  .miry+ 

^csm{~r 


m  ,miry+ .  .  ,nirz+  ' 
-ii,  j  cos(  — )  sin(  — ) 


n7T2T 

)cos(— — ) 


(A.10) 


(A. 11) 


Remaining  quantities  in  (A. 2)  and  (A. 3)  are 


7™  =  -  t2  (A.  12) 

and 

Z™  =  j^e  (A.  13) 

In  (A. 12),  k  —  uiy/jle  is  assumed  to  be  a  real  wave  number.  Here,  u  is 
the  angular  frequency  and  fi  and  e  are,  respectively,  the  permeability  and 
permittivity  of  the  homogeneous  medium  inside  the  waveguide.  The  radicand 
in  (A.  12)  is  therefore  purely  real  so  that  7 mn  is  either  purely  real  or  purely 
imaginary.  If  7mn  is  purely  real,  we  take  7mn  >  0.  If  7 mn  is  purely  imaginary, 
we  take  the  imaginary  part  of  7mn  to  be  non-negative. 
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From  [4,  sec.  8-1],  we  have 


=  £f(s-+,^+)e-’”-1 

f£?+  =  Y™&S(y+,z+)e 


"7mnX 


+  14. 


fcmnV’mf(y  +  ,'J+)e  1,m"X 


and 


Hi*'  =  -lSfiS(»+.*+)«~  +  u. 


sTElTE 


JUJU 


k2mn'Pin(y+,Z+)e'ymnX 

jun 


where 


4!2(»+.*+)  =  -»©si+,*+) 

The  wave  function  V>mf  (y+> 2+)  satisfies 

vvs; + = o 

subject  to  the  boundary  condition 

b.-E*”-0 

on  the  walls  of  the  waveguide.  Here,  un  is  the  unit  vector  normal 
of  the  waveguide.  Solutions  to  (A. 18)  and  (A. 19)  are 


mn 

V&£(v+,*+)  = 


^n  =  (^)2  +  (-)2 

0  C 


1  l^m^n 

be 


■ 


cos( 


■)  cos( 


nxz 


m  =  0, 1, 2 

>  n  =  0,1,2. 

m  +  n/  0 


where  en  is  Neumann’s  number  given  by 

f  1,  n  =  0 
"  \  2,  n  =  1,2,--- 

The  preceding  is  normalized  so  that 

Jo  dy+  Jo  dz+  {^^"(^+’z+)}  '  {2V’mf(y+^+)}  =  1 


j  (A. 14) 

■  (A. 15) 

A 

(A. 16) 
( A.17) 

(A. 18) 

(A. 19) 

to  the  wall 

•••  (A. 20) 

(A. 21) 

(A. 22) 
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Substituting  ip™  °f  (A. 20)  into  (A. 16)  and  (A. 17)  and  taking  2  with  respect 
to  the  coordinates  y+  and  z+ ,  we  obtain 


TEt  +  +n  T  (tmGi  f  n  ,mxy+x  .  ,nirz+ 
'*  >  "  («»-cos(— )sm(  — 


m  .  ,miry+ .  nirz+.' 
“It.  J  sm(  — ~ )  cos(  — — ) 


lTE/  +  +x  *  /^Tf  m  •  ,mxy+,  ,n5T2  + 

^  ' 2  >  =  kZj—  t“»T sm<— > cos(— 


n  .miry+ .  .  ,nirz+ . 

+U2-  cos(  — - — )  sin( - ) 

CO  c 


(A. 23) 


(A. 24) 


Remaining  quantities  in  (A. 14)  and  (A. 15)  are  7mn  given  by  (A. 12)  and 


yTE 
*  mn 


7mn 


(A. 25) 


Orthogonality  relationships,  are 

f  [  **(*!•.  ■  4,)  =  jf  jf  ^+(^.n  ■  &) 

=  /*<(»+  x4si-u  =1 11  ,Am’n)  =  (4- 

7o  y  io  ***'  “*  \  0,  otherwise 


(A. 26) 


where  0  is  either  TM  or  TE  and  6  is  either  TM  or  TE.  Equality  of  the 
three  integrals  in  (A. 26)  follows  from  (A. 4),  (A. 5),  (A. 16),  and  (A. 17).  It  is 
evident  from  [4,  eq.  (8-37)],  (A. 4),  and  (A. 9)  that  the  integrals  are  equal  to 
the  right-hand  side  of  (A. 26)  when  0  =  6  =  TM.  It  is  evident  from  [4,  eq. 
(8-36)],  (A. 17),  and  (A. 22)  that  the  integrals  are  equal  to  the  right-hand  side 
of  (A. 26)  when  0  =  6  =  TE.  Because  of  [4,  eq.  (8-38)],  the  integrals  vanish 
when  0^8.  If  there  is  a  degeneracy,  that  is,  if  kmn  =  kpq  for  (m,  rt)  (p.  q), 
then,  as  pointed  out  in  [4,  p.  390],  we  must  have 

t  iy*  f  dz*  Ut,(y* =  0,  P  =  TE,  TM  (A.27) 

JO  JO 

in  order  for  (A. 26)  to  hold.  It  can  be  shown  that  ip™  of  (A. 8)  and  ip™  of 
(A. 20)  satisfy  (A.27). 
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In  Appendix  A,  we  found  the  four  kinds  of  modes  that  can  exist  in  the 
rectangular  waveguide.  There  are  two  kinds  of  TM  modes  and  two  kinds  of 
TE  modes.  The  TM  modes  are  (E™+ ,  )  of  ( A.2)  and  (E™~ ,  ) 

of  (A.3).  In  (A.2)  and  (A.3),  s™,  and  hSxn  ***  Siven  by  (A-8),  (A. 10), 

and  (A. 11),  respectively.  The  TE  modes  are  (££f+,i££f+)  of  (A. 14)  and 
(£S?",2Z£?~)  of  (A. 15).  In  (A. 14)  and  (A.15),  and  are  8iven 

by  (A.20),  (A. 23),  and  (A. 24),  respectively. 
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Appendix  B 

Modes  of  the  Circular 
Waveguide 


Eight  kinds  of  modes  can  exist  in  the  circular  waveguide  whose  cross  sec- 

?TMe+  trTMe+ \  j  pTMe-  uTMe—  \ 
■np  >  ^-np 
?TEe 


tion  is  shown  in  Fig.  B.l.  They  are  (f£™e\iC[pMe+),  {E^e~ ^e~), 


/  jpTMo+  tjTMo+  \  /  rpTMo—  tjTMo—  \  /  r?TEe+  rjTEe+\  /  r?TEe-  rjTEe-\ 

\^np  -np  ),  \£Lnp  i  *Lnp  V^np  >  ^np  )i  V^np  ’  ±*-np 

{e£v°+,K!£0+),  and  Here,  R  is  the  electric  field  and  K 

is  the  magnetic  field.  The  superscript  “TM”  denotes  transverse  magnetic, 
uTEn  denotes  transverse  electric,  “e”  means  even  in  0,  “o”  means  odd  in  <f>, 
“+”  indicates  that  the  wave  travels  in  the  +z  direction,  and  ”  indicates 
that  the  wave  travels  in  the  —z  direction.  Here,  <j)  is  the  angle  measured 
counterclockwise  from  the  positive  x  axis  and  z  is  the  cylindrical  coordinate 
measured  in  the  up  x  direction  where  and  are  the  unit  vectors  in 
the  p  and  <f>  directions,  respectively.  By  definition,  p  =  \Jx2  4-  y2. 

From  the  analysis  in  [4,  sec.  8-1],  we  obtain 


Kv 


TAfe+  _  yTMeo  TMe/ 


ryl  ivj  CO  l  M«/  ^  IN. 
Znp  &np  (P,0)e 


-1™Z 


+  li* 


(fc™)2C(^)e‘ 


nip 


TMe+  _  lTMc 


JM, 


JUt 


(B.l) 


and 


TM  , 


K 


TMc- 


=  -z, 


TMeo  -TMc 


TM, 


TMc- 


np 
TMet 


^np 


■-(p,<j>)e^*  +  uz 


(k™)2ij;™'(p,  4>)e^p 


-np 


=  MpMe(p, 


JUt 


(B.2) 
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Figure  B.l:  Cross  section  of  the  circular  waveguide. 


where  u,  is  the  unit  vector  in  the  z  direction.  Moreover, 

C'(P'*)=  (B.3) 

«  =  -a,  X  YV>™‘(/>, *)  (B.4) 

The  wave  function  is  a  0)  which  is  even  in  4>  and  which 

satisfies 

VV™(/>,  #)  +  (*™)V™(/>, *)  =  0  (B.5) 

subject  to  the  boundary  condition 

0™  =  0,  p  =  o  (B.6) 

Even  solutions  to  (B.5)  and  (B.6)  are 


n  =  0,1,2,  ••• 
p  =  1,2,3,-  •  • 


(B.7) 
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where  Jn  is  the  Bessel  function  of  the  first  kind  of  order  n  and  {0  <  xnl  < 
xn2  <  Zn3  •  •  •}  satisfies 


Jn(xnp )  —  0,  p  —  1, 2,  3,  •  ■  * 

Furthermore,  e„  is  Neumann’s  number: 

f  1,  n  =  0 
n  \  2,  n  =  1,2,  •  •  • 

The  preceding  V*™*  is  normalized  so  that 

j r pdp  £’ ■  {s€p"ete«}  =  l 

To  verify  (B.10),  we  first  use  (B.7)  to  obtain 

so*  *)  =  cos(n« 


(B.8) 


(B.9) 


(B.10) 


sin  (n<j>) 
*&"/> 


}  (B-ll) 


If  the  left-hand  side  of  (B.10)  is  called  7lt  then  substitution  of  (B.ll)  gives 


ica2JZ+l(xnp) 


J  pdp  J  d<f> 


f  r/2  ,  ^n(CV)  sin2  {ruf>)  } 

■  (*nP  /»)C0S  (n<j>)  +  - (k™Pj2 - J 

Evaluation  of  the  integral  with  respect  to  <f>  reduces  (B.12)  to 


a2J2+1(Inp) 

Substituting 

into  (B.13),  we  obtain 


(B.12) 


I  J'2(k™p)  +  n2J^k™ 

l " 1  "p  p)  (k™p)2  J 

(B.13) 

x  =  k™p 

(B.  14) 

2/2 

(B.15) 

XnpJn+liXnp) 
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where 

h  =  Jq  '  y:J'n2(x)  +  ^Jl(x)Jdx 

It  is  shown  in  Appendix  C  that 

jf  {l-CV)  +  dx  = 

+jJ'„2(d)  +  dJn(d)J'„(d) 

Replacing  d  by  xnp  in  (B.17)  and  using  (B.8),  we  obtain 


Now  [6,  formula  9.1.27], 

J'n(x)  =  ~Jn+i(x)  +  -Jn(x) 

X 

Replacing  x  by  xnp  in  (B.19)  and  using  (B.8),  we  obtain 

A(Xnp)  =  ~Jn+  l(Xnp) 


(B.16) 


(B.17) 


(B.18) 


(B.19) 


(B.20) 


so  that  (B.18)  becomes 


(B.21) 


Substitution  of  (B.21)  into  (B.15)  gives  I\  =  1.  Thus,  (B.10)  is  verified. 
Substitution  of  (B.ll)  into  (B.3)  and  (B.4)  gives 


TMt 

~np 


(B.22) 


(B.23) 
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Remaining  quantities  in  (B.l)  and  (B.2)  are 


(B.24) 


and 


~TM 

yTMto  _  l np 

np  ~  jut- 


(B.25) 


In  (B.24),  k  =  Uy/jH  is  assumed  to  be  a  real  wave  number.  Here,  ui  is 
the  angular  frequency  and  p  and  e  are,  respectively,  the  permeability  and 
the  permittivity  of  the  homogeneous  medium  inside  the  waveguide.  The 
radicand  in  (B.24)  is  therefore  purely  real  so  that  7™  is  either  purely  real 
or  purely  imaginary.  If  7™  is  purely  real,  we  take  7™  >0.  If  7™  is 
purely  imaginary,  we  take  the  imaginary  part  of  7™  to  be  non-negative. 
The  superscript  “eo”  was  placed  on  the  left-hand  side  of  (B.25)  to  avoid 
confusion  with  the  quantity  Z™  which  was  defined  in  Appendix  A. 

Similar  to  (B.1)-(B.4),  we  have 


r?TMo+ 

^np 


TMo+ 


_  yTMeo  TMo 

^ np  “np 


(p,  +  uz 


JUS. 


=  4i 


TMo 

np 


jut 


(B.26) 


and 


where 


Y  Jl?°(pA)  (B.28) 

-!1.xVO’»)  (B.29) 

The  wave  function  0™°(p,<£)  is  a  ij>™{p,  <t>)  which  is  odd  in  o  and  which 
satisfies  (B.5)  and  (B.6).  Odd  solutions  of  (B.5)  and  (B.6)  are 


TMo 
^np 


uTMo 

— np 


=  - 
(p,<f>)  =  - 


0)e7™J 


+  MZ- 


KM)2*™°(p,<t>)e 


,TM  , 


np 


JUt 


(B.27) 


k™ 

np 


K  (P<  *) 


F2Jn(k™  p)sin(n<t>) 

7T 


np 

xnpi/, 


n  =  1,2,3, 
P=  1,2,3, 


n+1 


(Xnp) 


(B.30) 
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The  preceding  0™°  is  normalized  so  that 


£pdp  ■  {EiC \pa)}  =  i  (b.31) 


To  verify  (B.31),  we  first  use  (B.30)  to  obtain 


sC °(p 


(Xnp) 


\  7#/»TM  x  .  /  IX  .  n^n(^p)cos(n0)] 

|Mn(*np  />)  Sin(n^)  +  2i4, - P-jjJTp - 1 


(B.32'1 


Substituting  (B.32)  for  we  find  that  the  left-hand  side  of  (B.31)  is 

the  right-hand  side  of  (B.12)  with  e„  replaced  by  2  and  with  sin(nd>)  and 
cos(n<£)  interchanged.  Integrating  this  result  with  respect  to  <f>,  we  discover 
that  the  left-hand  side  of  (B.31)  is  Ix  of  (B.13).  Now,  as  stated  in  the  sentence 
following  (B.21),  I\  =  1.  Thus,  (B.31)  is  verified. 

Substitution  of  (B.32)  into  (B.28)  and  (B.29)  gives 


^  V*  \aJn+l(xnp) 


)  rl/lTM  \  ■  /  n  ,  n^n{knp 

{UpJniKp  p)sm{n4>)  +  u^ - £ 


nJn{k™ p)  cos(n<t>) 
k™P 


(B.33) 


Kp  °(P,  <t>)  V  7T  U«/n+l(xnp) 


f  n Mk™p)  cos{n<t>)  TlftTM  \  ,00<N 

|Mp - - u*Jn(knp  P)  sin(n<p)|  (B.34) 


The  quantities  and  Z™eo  jn  (B.26)  and  (B.27)  are  given  by  (B.24)  and 
(B.25),  respectively. 

From  [4,  sec.  (8-1)],  we  have 


KTJ '*  =  YLE‘°h™'(p,  *)e-*'*  +  u: 


(knp)2^npe(P^)e  ^npZ  f  (B-35) 
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and 


where 


£p'(p,4>)  ~  U.z  X  Yrp™e(p,<t>) 

tf?W)  =  -vC(p^) 


The  wave  function  tpnp'iPi  <j>)  is  a  ip[p(p,<f>)  which  is  even  in  <j>  and 

QJltlCn  po 

vV2>,  *)  +  (*Sf )  V«  (*  *)  =  0 

subject  to  the  boundary  condition 


dp 


=  0, 


Even  solutions  of  (B.39)  and  (B.40)  are 


uTE  _  fza 
_  ■“** 

np  a 


MknpP)cos{n<f>) 


.  ^'np  -  ”2)  */««p) 

where  {0  <  x(,x  <  x^  <  x^  ■  •  •}  satisfies 


n  =  0, 1, 2,  •  •  • 
p  =  1,2,3,  ••• 


J'n«p)  =  0,  p  =  1,2,3,--* 

The  preceding  is  normalized  so  that 

[pdf'  [’  dd>  {Z4>1?(p,  «}  •  {2iC'(*  *)} 

To  verify  (B.43),  we  first  use  (B.41)  to  obtain 


=  1 


(B.36) 

(B.37) 

(B.38) 

which 

(B.39) 

(B.40) 

(B.41) 

(B.42) 

(B.43) 

(B.44) 
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If  the  left-hand  side  of  (B.43)  is  called  /3,  then  substitution  of  (B.44)  gives 


h  = 


L.TE 

K"P 


*{xnl~n2)  \MXnp) 


(■ 


J'n2(klpP)  cos 2{n<t>)  + 


£pdpjQ  d<t> 
^Jlifipp)  sin  2{n<f>)' 


(k™p)2 


(B.45) 


Evaluation  of  the  integral  with  respect  to  <f>  reduces  (B.45)  to 


h  = 


k™ 


<1  -  «2  \M<P) 

Substituting 


into  (B.46),  we  obtain 


h  = 


x  =  klpP 
2A 


where 


u 


=  £"{xfn\x)  +  n2JMyx 

Replacing  d  by  x'np  in  (B.17)  and  using  (B.42),  we  obtain 

r  «2  -  n2)72(x,np) 

U  -  n 


(B.47) 


(B.48) 


(B.49) 


(B.50) 


Substitution  of  (B.50)  into  (B.48)  gives  /3  =  1.  Thus,  (B.43)  is  verified. 
Substitution  of  (B.44)  into  (B.37)  and  (B.38)  gives 


<p^(p,4>)  = 


L.TE 

Knp 


7r(x;2  _n2)  \Jn(x'np) 
njn{k™p)  sin (n<£) 


{a- 


*JF/> 


+  U^Jnik^p  p)  cos{n<t>)^  (B.5i; 


K?M)  =  - 


lTE 

—HE— 


{ 


*(x'np  ~  "2)  \Mxnp)J 

ti  1 1 TE  \  /  njn{k™p)sm{n0)} 

UpJM; '  p)  cos{n<j>) - 1  (B.o2) 
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Remaining  quantities  in  (B.35)  and  (B.36)  are 


iTJ  =  MFy  -  * 


(B.53) 


yrTEeo  __  >np 

np  jun 


(B.54) 


The  radicand  in  (B.53)  is  purely  real  so  that  7^  is  either  purely  real  or  purely 
imaginary.  If  7^  is  purely  real,  we  take  7^  >  0.  If  7^  is  purely  imaginary, 
we  take  the  imaginary  part  of  7^  to  be  non-negative.  The  superscript  “eo” 
was  placed  on  the  left-hand  side  of  (B.54)  to  avoid  confusion  with  the  quantity 
Y™  which  was  defined  in  Appendix  A. 

Similar  to  (B.35)-(B.38),  we  have 


(kIv)2^np°(P^)e  ~'npZ  f  (B‘55) 


2  f  (B.56) 


where 


’(**)  =  a.*  2*  -TfW) 

’(/>,  <f>)  =  0) 


(B.57) 

(B.58) 


The  wave  function  tp^°(p,<f))  is  a  4^p{p,4>)  which  is  odd  in  o  and  which 
satisfies  (B.39)  and  (B.40).  Odd  solutions  of  (B.39)  and  (B.40)  are 


LTE  _  nP 

“  a 


0 lp°(P » <A)  = 


2  4(^np£/?)siD(n<^) 

*«p-"2)  ^«p) 


n  =  1,2,3, 
P=  1,2,3.-  •• 


(B.59) 
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(B.60) 


The  preceding  ipnp° 's  normalized  so  that 

J0  PdpjQ  d<f>  {WnP£°(P>  <£)}  •  {Y rl>™°{p,  <£)}  =  1 

To  verify  (B.60),  we  first  use  (B.59)  to  obtain 


Substituting  (B.61)  for  Y^f0 ,  we  find  that  the  left-hand  side  of  (B.60)  is 
the  right-hand  side  of  (B.45)  with  tn  replaced  by  2  and  with  sin (n<f>)  and 
cos(n<f>)  interchanged.  Integrating  this  result  with  respect  to  <f>,  we  discover 
that  the  left-hand  side  of  (B.60)  is  /3  of  (B.46).  Now,  as  stated  in  the  sentence 
following  (B.50),  I3  =  1.  Thus,  (B.60)  is  verified. 

Substitution  of  (B.61)  into  (B.57)  and  (B.58)  gives 


sin(n0)  +  34- 


nMknpp)  cos{n<t>) 


kTnpEP 


(B.62) 


(B.63) 


The  quantities  7^  and  YjpEeo  in  (B.55)  and  (B.56)  are  given  by  (B.53)  and 
(B.54),  respectively. 

Orthogonality  relationships  are 


£'«<*•  *5)  .. 


1,  (t,r,m,q)  =  {u,s,n,p) 
0,  otherwise 


;B.64) 
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where 


r  =  e,  o 
s  =  e,  o 
t  =  TE,  TM 
u  =  TE ,  TM 

Equality  of  the  three  integrals  in  (B.64)  follows  from  (B.3),  (B.4),  (B.2S), 
(B.29),  (B.37),  (B.38),  (B.57),  and  (B.58).  It  is  evident  from  [4,  eq.  (8-37)], 
(B.3),  (B.10),  (B.28),  and  (B.31)  that  the  integrals  are  equal  to  the  right- 
hand  side  of  (B.64)  when  t  =  u  =  TM  and  r  =  s.  When  t  =  u  =  TM 
and  r  s,  the  integrals  vanish  because  the  integrands  are  odd  in  <p.  It  is 
evident  from  from  (4,  eq.  (8-36)],  (B.37),  (B.43),  (B.57),  and  (B.60)  that 
the  integrals  are  equal  to  the  right-hand  side  of  (B.64)  when  t  =  u  =  TE 
and  r  =  3.  When  t  =  u  =  TE  and  r  s,  the  integrands  vanish  because 
the  integrands  are  odd  in  <f>.  Because  of  [4,  eq.  (8-38)],  the  integrals  vanish 
when  t  ^  u. 

In  Appendix  B,  we  found  the  eight  kinds  of  modes  that  can  exist  in  the 
circular  waveguide.  There  are  two  kinds  of  even  TM  modes,  two  kinds  of 
odd  TM  modes,  two  kinds  of  even  TE  modes,  and  two  kinds  of  odd  TE 
modes.  The  two  kinds  of  even  TM  modes  are  ,  H™e+ )  of  (B.l)  and 

(£["*-, of  (B.2).  In  (B.l)  and  (B.2),  v>™',  si"',  and  /£"'  are 
given  by  (B.7),  (B.22),  and  (B.23),  respectively.  The  two  kinds  of  odd  TM 
modes  are  (£™^£™°+)  of  (B.26)  and  (S™-,#™-)  Gf  (B.27).  In 
(B.26),  and  (B.27),  e^°,  and  are  given  by  (B.30),  (B.33),  and 

(B.34),  respectively.  The  two  kinds  of  even  TE  modes  are  (E^pe+  ,Knpe+)  of 
(B.35)  and  )  of  (B.36).  In  (B.35)  and  (B.36),  <£fe,  and 

are  given  by  (B.41),  (B.51),  and  (B.52),  respectively.  The  two  kinds  of 
odd  TE  modes  are  (i£p£o+, l£pEo+)  of  (B.55)  and  (^°~ , KTnpEo~ )  of  (B.56). 
In  (B.55),  and  (B.56),  0^pEo,  e^°,  and  h™°  are  given  by  (B.59),  (B.62),  and 
(B.63),  respectively. 
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Appendix  C 

Evaluation  of  an  Integral  of 
Bessel  Functions 


In  Appendix  C,  we  evaluate  the  integral  I  defined  by 


I  =  £  {<’(*)  +  dx 

(C.l) 

Proceeding  as  in 

[7,  sec.  5.296],  we  recast  (C.l)  as 

(C.2) 

Substitution  of  [6,  formula  9.1.27] 

2n 

*/n(®)  =  wfn_i(x)  -f"  </n-(-i(x) 

X 

(C.3) 

and 

2 J'n{x)  =  Jn-i{x)  -  Jn+ i(x) 

(C.4) 

into  (C.2)  yields 

/  =  +  /„+I) 

(C.5) 

where 

/n  =  /  lJl{x)dx 

Jo 

(C.6) 
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To  evaluate  /„  of  (C.6),  we  will  use  Bessel’s  equation  [6,  formula  9.1.1]: 


(C.7) 


(C.8) 


(C.9) 


x2J"(x)  +  x  J'n ( x )  +  (x2  -  n2)  Jn(x)  =  0 
Multiplying  (C.7)  by  ./'(x),  we  obtain 

i >JW'(x)  +  xj'„\x)  +  (iJ  -  »V.(*MU*)  =  0 

Expression  (C.8)  is  recast  a s 

j  (42M)'  +  +  j  (4M)'  -  j  (4M)'  =  o 

The  integral  of  (C.9)  from  0  to  d  with  respect  to  x  is 

5  Jo  **  (4J(*))'  *  +  jo  zJ,»‘(z')dx  +  \Ja  z2(j2(z))'dz 

-y  (4(<0  -  4(0))  -  o  (C.io) 

Because  [6,  formula  9.1.7] 

-^(0)  =  {  „  =  1,2,3, 


(C.ll) 


the  J%(0)  term  drops  out  of  (C.10).  Integrating  the  first  and  third  integrals 
in  (C.10)  by  parts,  we  obtain 


j  (4!(<<)  +  4(<0)  -  £ *4(*)<k -  y4(<0  =  o  (C.i2) 

whence 

/n  =  \  +  (d1  -  n!)4(<0}  (0.13) 

Substitution  of  (C.13)  into  (C.5)  gives 

/  =  i{(«f  -  (n  -  1  )24_,(<1)  +  (J1  -  (n  +  l)J)4+iM 

+J2Al,(d)  +  dXU(d)}  (C.14) 
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Adding  (C.4)  to  (C.3),  we  obtain 

Jn-i{x)  =  —Jn(x)  +  J'n{x)  (C.15) 

x 

Subtracting  (C.4)  from  (C.3),  we  obtain 

Jn+1(x)  =  -Jn(x)  -  J'n(x)  (C.16) 

x 

From  [6,  formula  9.1.27],  we  have 

JL  ,(*)  =  -W  +  (C.i7) 

X 

Jn+l(x)  =  Jn{x) - Jn+\(x)  (C.18) 

X 

Substitution  of  (C.15)  into  (C.17)  gives 

•£-.(*)  =  (“I  +  ■«*)  +  (C.19) 

Substitution  of  (C.16)  into  (C.18)  gives 

•/;+.(*)  =  (l  -  — *(*)  +  (C.20) 


Using  (C.15),  (C.16),  (C.19),  and  (C.20)  to  convert  the  Bessel  functions  of 
order  n  ±  1  in  (C.14)  to  Bessel  functions  of  order  n,  we  arrive  at 

I  =  J< S  -  n2)Jl(d)  +  jW)  +  dJ,(i)J‘M)  (C.21) 

This  result  agrees  with  [7,  eq.  (12)  on  page  186]. 
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Appendix  D 

Radiation  of  a  Magnetic 
Current  Element  in  a  Circular 
Waveguide  Matched  at  Both 
Ends 


In  Appendix  D,  we  determine  the  fields  £(Q, —  n')),  i£(Q,  146(1:  —  z:')), 
£(Q,  U.zS(l  —  l')),  and  if(Q, —  n'))  in  the  circular  waveguide  whose  cross 
section  is  shown  in  Fig.  B.l.  For  simplicity,  we  assume  that  the  waveguide  is 
matched  at  both  ends.  Here,  K{d.v,M.v )  is  the  electric  field  due  to  the  combi¬ 
nation  of  the  electric  current  source  J?  and  the  magnetic  current  source  Mv 
where  both  J?  and  M 7  are  volume  densities  of  current.  Similarly,  K{JLV,ALV) 
is  the  magnetic  field  due  to  (iZv,Mv)-  The  argument  “Q”  indicates  that  there 
is  no  electric  current  source  in  the  waveguide.  Furthermore,  6(r  —  r')  is  the 
three  dimensional  Dirac  delta  function,  is  the  unit  vector  in  the  cj>  direc¬ 
tion,  and  u.z  is  the  unit  vector  in  the  z  direction.  Here,  r  is  the  radius  vector 
to  the  point  whose  cylindrical  coordinates  are  ( p,<j>,z ),  and  r'  is  the  radius 
vector  to  the  point  whose  cylindrical  coordinates  are  {p\  <j>' ,z'). 

Consider  £.(Q,HpS(r  —  r'))  and  i£(Q,  !4<5(r  —  r'))  evaluated  at  the  point 
whose  cylindrical  coordinates  are  (p,  <f>,  z).  Since  the  only  source  is  a  trans¬ 
verse  magnetic  current  located  at  (/?',  <{>',  z')  and  since  the  waveguide  is  matched 
at  both  ends,  there  will  only  be  waves  that  travel  outward  from  z'  and  the 
transverse  component  of  the  magnetic  field  will  be  continuous  at  z' .  Hence, 
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£(Q,24^(£  ~  L1))  and  2£(Q,  246(1:  —  l'))  may  be  expanded  as 

(fitWte-E'W.fitWfe -£')))  =  E  E  E 


n=0  p= 1  *=e,o 


±  t,£Sfa±)}  (D.l) 

where  the  upper  sign  is  to  be  taken  for  z  >  z'  and  the  lower  one  for  z  <  z'. 
In  (D.l),  (£™*±,C",±)  and  (£^f ** )  are  given  by  (B.l),  (B.2), 
(B.26),  (B.27),  (B.35),  (B.36),  (B.55),  and  (B.56).  Moreover,  7™  and  7jp£ 
axe  given  by  (B.24)  and  (B.53),  respectively.  The  s  =  o  term  is  to  be  omit¬ 
ted  from  (D.l)  when  n,  the  index  of  the  outer  summation,  is  zero.  In  (D.l), 
CTMa,M<t>  ajjd  (jTEa,M<t>  ajg  unknown  constants.  The  superscript  UM<£”  indi¬ 
cates  that  the  source  is  the  <f>  directed  magnetic  current  246(2;— r').  Substitut¬ 
ing  the  expressions  in  Appendix  B  for  (E™a± ,  H™3± )  and  ,  R^3± ) 

into  (D.l),  we  obtain 

£(Wte-n'))  =  EE  E  c™"Mt 

n=0  p=l  *—e,o 


(e'wwv 


T»=0  P=1  a=e,0 

m ,!!*%-/))  =  EE  E 

n=0  p=l  J=«,o 

+  EE  E  ci£''"*{^E”4n,E'(/>,*) 

n=0  p=l  J=e,o 

e(^,^)(^pE)20ffJ(p^)  | c-,r« 

jun  J 


(D.2) 


where 


Now, 


.  f  1,  z  >  z' 

'<*’*>  =  |  -1,  z<z' 


S(l  -  r')  =  S(£  -  £)^(2  -  2') 


(D.3) 


(D-4) 


(D.5) 
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where  £  is  the  radius  vector  to  the  point  in  the  xy  plane  whose  cylindrical 
coordinates  are  (p,  4>)  and  £'  is  the  radius  vector  to  the  point  in  the  xy  plane 
whose  cylindrical  coordinates  are  (p1,  <j>').  Furthermore,  6(p  —  £')  is  the  two- 
dimensional  Dirac  delta  function,  and  S(z  —  z')  is  the  one-dimensional  Dirac 
delta  function.  Thanks  to  (D.5),  the  fields  on  the  left-hand  sides  of  (D.2)  and 
(D.3)  become  £(0.  M8(z  —  z'))  and  H( 0.  M8(z  —  z1))  where  M  is  a  surface 
density  of  magnetic  current  in  the  z  —  z'  plane: 

M  =  -£')  (D.6) 

The  right-hand  side  of  (D.2)  is  discontinuous  at  z  —  z' .  The  discontinuity  is 
related  to  M.  by  [4,  eq.(3.14)] 

K  =  (R+  -  R~)  x  Hz  (D.7) 

where 

£+  =  Km£(Q,Atf(*-*')) 

z—z' 
z>z' 

ET  =  lim 

z-*z' 
z<z' 

Taking  u*x  of  both  sides  of  (D.7),  we  obtain 

u,  xM  =  «z  x  ((£+-£')  xuz)  (D.10) 

or,  more  simply, 

U.z  *  M.  =  {£L+  ~  E.  )tan 

where  the  subscript  “tan”  denotes  the  transverse  component 
(D.ll)  to  the  discontinuous  electric  field  of  (D.2),  we  obtain 

-M(2-2')  =  2 EE  E 

n=0  p=l  s=e,0 

+2  E  E  E  *)  (D12) 


(D.ll) 

Applying 


(D.8) 

(D.9) 
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Scalar  multiplication  of  (D.12)  by  §*q{p,  <t>)  and  integration  o^er  the  cross 
section  of  the  waveguide  give 

-Jo  pipL  ^  W s(£- e') 

=  2  E  E  E  f  pdp  r  H  {^,(/,,  *>} 

n=0  p=l  *=e,o  J0 

+2±±  E  C™‘M*  [pip  /2'<W{^((.,«  •*"*((., «}(D.13) 

n=0  p=l  i=e,o  -70  1/0 


In  (D.13),  we  choose 


r  =  e,o 


t  =  TM,TE 
m  =  0,1,2,  ••• 
9  =  1,2,3,--- 


(D.14) 


The  definition  of  £(£—£')  and  the  orthogonality  relationships  (B.64)  are  used 
to  evaluate  the  integrals  in  (D.13).  Next,  r,  m,  and  q  axe  replaced  by  s,  n, 
and  p ,  respectively.  The  result  is 


^TMs,M4>  V  '  £np  *(/A 

^nP  2  Z™eo 

where  up,  is  the  unit  vector  in  the  /?'  direction. 

Substitution  of  (D.15)  and  (D.16)  into  (D.2)  and  (D.3)  gives 

m, iuKt  -  r'))  =  1 E  E  (v '  4"V,  *')) 

"  n=0  p=l  4=e,o 

f./_  JTMsf  _  ^  ,  ..  (*npM)2^pM,(/?,^)l 


(D.15) 

(D.16) 


,/  _  _/\  JTMsf  _  ± \  ,  _  v  np  7  ^rip  W 

«(*’*)**  - ]UeZ™'° 


e-~lnp  \z~z 


-sEE  E  <D17) 


n=0  p=l  4=e,o 


£(Q> 14%  -  cO)  =  -5  £  £  H 

“  n=0  p=l  *=e,c 


(v  •  jj) 


7T  Mco 

np 
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-5EEE(v  Cv,  *')) 

n=0  p=l  3=e,o 

•  (d.18) 

Substituting  (B.22),  (B.33),  (B.51),  (B.62),  (B.7),  (B.30),  and  (B.25)  into 
(D.17),  we  obtain 


R (a, -r'))  =  iEE -f  -  2  7  np, 

^  n=Op=l  Inp‘^n+l(Inp) 


^(*5?)W™V)/ 
■—zrjrTn — l 


■J'n{klp  p)  cos(n(<t>  -  0'))  +  1 ^ 


e(2,  z')nJn(k™p )  sin(n(<p  -  o')) 


1  "  "  g(gt^)(fcff)»nJn(fc^) 
*  «2  ~  n2)(A£pV)^(x'„p) 


,  u  fcnpM-/n(^pMp)  cos(n(0  -  ^))\ 

-yTAf  }e 

I np 

1  ”  ”  t(z,z')(klr)2nJn(k™p')  /  nJJk^p) 

»  S  &  «2p  ~  n2)(^pV)J2(x'np)  ^ 

•  cos(n(0  -  0'))  +  UjJ'niknp  p)  sin(n(0  -  0O)}e“7"p£|*'~2,1  (D.19) 

Similarly,  substitution  of  (B.22),  (B.23),  (B.33),  (B.34),  (B.51),  (B.52),  (B.62). 
(B.63),  (B.41),  ^B.59),  and  (B.54)  into  (D.18)  gives 

R( Q,^(r  — t)j  - -3—  EE  jmTjS  1> 

'57rU;A»  n=0p=l  'np  XnpJn+l>,xnpJ 

nj  (fc™p)  1 

- — sin(n(0  -  <?'))  +I4-/n(^npWp)cos(n(0  -  0'))} 


_ TM  |  _  i\ 

.g — I'np  I*-*  I 


;  (QW.^V)  r 


•7 ™A(k™p)  sin(n(4,  -  6'))  +  u, 


z'lklf  JJklf  p)  smiXe  -  ^'))}e 


*J?P 


T&fl*-*'! 


(D.20) 


Consider  R(Q.,u.z6(r  —  r'))  and  H_(Q.,iLz6{l  —  l'))  evaluated  at  the  point 
whose  cylindrical  coordinates  are  ( p,®,z ).  To  obtain  expressions  for  these 
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fields,  we  replace  the  longitudinal  magnetic  current 


Kv  =  uz6(r.  ~  r!) 

(D.21) 

by  an  equivalent  transverse  electric  current  [8].  To  see  that  such  a  “re¬ 
placement”  is  possible,  recall  that  ( E( 0.  Mv ).  H( Q.  Mv))  satisfies  Maxwell’s 
equations: 

£  x  R{Q,M?)  =  -j^m,Mv)M.v  \  m  w 

Y  x  m,M.v)  =  jueK(b,M.v)  J  1  ; 

If 

Mv 

K  =  K(bM.v)  +  =z- 

(D.23) 

then  (D.22)  becomes 

2x£(Q,r)  =  -wi  1 

£x£  =  jutmMv)+lv  J 

>  (D.24) 

where 

Vxtf 

jufi 

(D.25) 

If  we  can  find  the  electromagnetic  field  {K{Q.,M.V),K)  that  appears  in  (D.24), 
then  we  will,  of  course,  have  2£(Q,MV)  and,  from  (D.23),  K(Q.,M.V)  will  be 
given  by 

m,M.v)=K~—  (D-26) 

JUfl 

From  (D.24),  the  electromagnetic  field  (f£(Q,  M.V),E.)  is  the  field  radiated 
by  the  electric  current  source  J?  so  that  we  may  write 


£(a,iT)  =  £U\Q) 

K  =  KUV,  Q) 


(D.27) 

(D.2S) 


Now,  our  objective  is  to  find  the  electromagnetic  field  (£(iZv,Q),i£(JZt\Q))- 
Unfortunately,  J_v  can  not  be  obtained  by  substituting  (D.21)  into  (D.25) 
because  S(r-r')  is  not  differentiable.  Substitution  of  (D.5)  into  (D.21)  gives 

Mv  =  uz6{z-z')8{p-£)  (D.29) 
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(D.30) 


We  attempt  to  define  A£u  by 


My  =  Uz6(z  -  z')  lim  /jv(£  -  p) 

N—oo  — 

where  /yv  becomes  more  and  more  impulsive  as  N  increases: 

=  S(E- (D.31) 

Now,  we  can  not  define  {E(Q.,MV),H.(Q.,MV))  by 

m,M.v)  =  -  zOJim  fN(&  -  £))  | 

m,M?)=  |  (D'32) 

N-+OO  -  j 


because,  since  £(£  —  is  not  a  legitimate  function,  the  limits  indicated  in 
(D.32)  do  not  really  exist.  The  proper  definition  of  (£(£),  A£u),  i£(Q,  Mv))  is 

E(Q,MV)  =  | 

=  -*')/*(£ -£))  /  (  33) 

Taking  the  approach  outlined  in  the  previous  paragraph,  we  replace  the 
fields  due  to  Mv  of  (D.21)  by  the  fields  due  to 


M.-  =  -.'lIEI  *)  (D.34) 

n=0  p=l  «=e,o 

where  the  s  =  0  term  is  to  be  omitted  when  n,  the  index  of  the  outer 
summation,  is  zero.  We  will  choose  C™a  such  that  Mv  of  (D.34)  approaches 
Mv  of  (D.29)  as  both  N\  and  N?  increase.  Finally,  we  will  take  lim,Vi,.v2— oo 
of  the  fields  due  to  A£v  of  (D.34). 

In  order  for  Mv  of  (D.34)  to  approach  M_v  of  (D.29),  the  triple  summation 
in  (D.34)  must  approach  6(p  —  Hence,  we  write 

c™’*zrE,(p, «)  (d.35) 

n=0  p=l  a =e,o 

Although  it  is  not  strictly  true,  (D.35)  is  a  means  of  evaluating  C^s .  Mul¬ 
tiplying  both  sides  of  (D.35)  by  T(P >  <P)t  integrating  over  the  cross  section 

of  the  waveguide,  using  the  orthogonality  relationship  [1,  secs.  8-1  and  8-2] 


(C  f  J°  P  ip  £'  i*  (p<  «)<?’(*>.  *) 


1,  (r,  m,  q)  =  (s,  n,  p) 

0,  otherwise 

(D.36) 
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(D.37) 


and  finally  replacing  (r,  m,q)  by  (s,n,p),  we  obtain 

cjf*  =  (*”  )V”V. « 

Substitution  of  (D.37)  into  (D.34)  gives 


(d-38> 

n=0  p=l  i=e,o 

Substituting  (D.38)  into  (D.25)  and  using  (B.38)  and  (B.58),  we  obtain 

lv  =  16(z  -  z')  (D.39) 


where 


ILt  x  ~ 


jun 


EE  E  (knp)2'l>Zp3(p\  <t>) 


n=0  p=l  »=e,o 


(D.40) 


Note  that  is  a  surface  density  of  electric  current  in  the  z  =  z'  plane. 

This  paragraph  and  the  next  paragraph  are  devoted  to  finding  the  electro¬ 
magnetic  field  {R(Q.,M.V),K)  which  satisfies  (D.24)  with  J?  given  by  (D.39). 
Since  the  only  source  is  a  transverse  electric  current  in  the  z  =  z'  plane,  and 
since  the  waveguide  is  matched  at  both  ends,  there  will  only  be  waves  that 
travel  outward  from  z'  and  the  transverse  component  of  the  electric  field  will 
be  continuous  at  z' .  Hence,  £(Q,MV)  and  H_  may  be  expanded  as 


(£(a,An,l)  =  EEE  {±c™‘Je* 


TM  „/ 


n=0  p=sl  s=c,o 

•(£!"’*.  alf'*)  +  *,i£ 


TEs± 

^~np 


where  the  upper  sign  is  to  be  taken  for  z  >  z'  and  the  lower  one  for  z  <  z' . 
In  (D.41),  ,££"■*)  and  (££,*•*, fljf**)  are  given  by  (B.l),  (B.2), 

(B.26),  (B.27),  (B.35),  (B.36),  (B.55),  and  (B.56).  Moreover,  7™  and  ' 
axe  given  by  (B.24)  and  (B.53),  respectively.  The  s  =  0  term  is  to  be  omitted 
from  (D.41)  when  n,  the  index  of  the  outer  summation,  is  zero.  In  (D.41), 
C™'"'  and  CLE‘J  are  unknown  constants.  The  superscript  J  indicates  that 
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the  source  is  ]_  of  (D.40).  Substituting  the  expressions  in  Appendix  B  for 
(£™3±,2tTM3±)  and  (RTEa±,RTEs±)  into  (D.41),  we  obtain 

=  ee  e 

J 


+  EE  E 

n=0  p=l  a=e,o 

■2  =  E  E  E 

n=0  p=l  *=e,o 

■e-’""1*-''  +  E  E  E  O'  {<(*,  *')0 


(D.42) 


(D-43) 


where  e(z,  z')  is  given  by  (D.4). 

The  right-hand  side  of  (D.43)  is  discontinuous  at  2  =  z'.  The  disconti¬ 
nuity  is  related  to  J_  of  (D.40)  by  [4,  eq.  (3-14)] 


1  =  U'*{K  -K) 


(D.44) 


where 


R  =  lim  H 


(D.45) 


R  =  lim  R 


(D.46) 


Taking  —  nzx  of  both  sides  of  (D.44),  we  obtain 


or,  more  simply, 


-  uz  x  J_=  -uz  x  ju,  x  ( R+  -  R  )  J 


-uz*l  =  (R  ~R)u 


(D.47) 


(D.48) 


where  the  subscript  “tan”  denotes  the  transverse  part  of  the  vector.  Applying 
(D.48)  to  the  discontinuous  magnetic  field  of  (D.43),  we  obtain 


-a.  x  i  =  2  E  E  E  4) 

n=0  p=l  J=e,o 

+2  E  E  E  *) 

rv=0  p=l  »=e,o 

Substitution  of  (D.40)  into  (D.49)  gives 

1  AT,  Nj 


(D.49) 


—  EEE  ( WC'V.  *')CW)  =  2 

JUr  n=Q  p-1  s=e,o 

•  E  E  E  ™-(p,  </>) + 2  f;  f;  £  <£**■ jy™‘°£*(P,  4)  (d.so) 

n=0  p=l  *=e,o  n=0  p=l  J=e,o 

Because  of  the  orthogonality  (B.64)  of  h^Ea  and  h™a,  it  is  evident  from 
(D.50)  that,  upon  using  (B.54)  to  dispose  of  YjpEeo, 

C™a'J  =  0  (D.51) 

f  (g)2g(/,^)  fSV^  v 

2-e  ’  1  (d.s2> 


rTEa,J  _ 
t-'np  ~ 


s  =  e,  o 

«  n  =  0, 1,  •  •  • ,  N\ 

,  p=  1,2, •••,Ar2 

otherwise 


(D.52) 


In  view  of  (B.54),  substitution  of  (D.51)  and  (D.52)  into  (D.42)  and  (D.43) 
gives 


m,m  =  i e E  E 

^  n=0  p=l  *=e,o  7np 

fl=-£lEE  ik™n™‘(p\4‘) 

n=0p=  1  J=e,o 


(D.53) 


f  /  /ur&z  IX  ,  ik™)2'l’np3(P,<P) 

<e(z,z)^  (p,  <j>)  +  uz — p  ~r| - 

l  'np 


(D.54) 


Substituting  (B.41),  (B.59),  (B.51),  and  (B.62)  into  (D.53),  taking 
lim/vj.jvj-oo,  and  then  substituting  (D.21)  for  M2 ,  we  obtain 

u  <5(r  -  r'))  =  -ff  Jn(knp)3Jn{knpp')  j  nJn(k™  p) 

£l{U,-zC{_  r  9  2^2^  TE(ti  2  _n2\T2(T>  \\^P  l.TE  n 


2x  „t-0ptl  7,rp£(x'4  -  n2)J2(x;p) 


*5fp 
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•sin (n(<j>  -  <f>'))  +  p)  cos(n((f>  -  <f>'))} e  ^z~z'\  (D.55) 

Substituting  (B.41),  (B.59),  (B.52),  and  (B.63)  into  (D.54),  we  obtain 

fj  i  ^  ^  tn(knf)3Jn(kZpP')  f  ,  ,s 


UaO^P)  cos( n(/p  -  4>'))  -  ji»— 


(klpp)  sin (n(<t>  - 

klfp 


klpMknpP)  cos(n(<£  -  <t>')) 


(D.56) 


We  want  to  substitute  K  of  (D.56)  and  M.v  of  (D.38)  into  (D.26)  before  letting 
Ni  and  Ni  approach  oo.  Substitution  of  (B.41)  and  (B.59)  into  (D.38)  gives 

M v_,  6(Z  -  Z')  ^  ^n(knf)2Mknp  P')Mknp  P)  COs{n{(f>  -  4>')) 

^ z  *  2-(  2-^  t-n  _  ni„i  \  (U.o t ) 


n=0  past 


«»  ~  n2)^«P) 


Substituting  (D.56)  and  (D.57)  into  the  right-hand  side  of  (D.26),  taking 
lim^jv^oo,  and  then  substituting  (D.21)  into  the  left-hand  side  of  (D.26), 
we  obtain 

r  „n  C,  rn  j  ^  ^  tn(knf)3Jn(knp  p')  /  /  ,  „ 


HpJ'n(klpp)  cos(n(d>  -  0'))  -  It 


nJ*(kZpP)  sin(n(<?>  -  <£')) 


knp'  P 


knpMklpP)  cos{n{<p  -  4>')) 


26(z  -  z’)Jn(knpP )  COS (n(<t>  -  O')) 


(D.58) 


There  axe  two  r±z  terms  inside  the  double  summation  on  the  right-hand 
side  of  (D.58).  The  second  uz  term  came  from  the  — Mv/(iun 1  term  on  the 
right-hand  side  of  (D.26).  Therefore,  if  we  deleted  the  second  uz  term  from 
the  right-hand  side  of  (D.58),  we  would  have  H_  instead  of  i£(Q,  uz8(r  —  r')). 
Equation  (8)  in  [9,  sec.  22]  correctly  gives  £(Q, U.z6(l  —  r'))  of  (D.55)  for  the 
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electric  field  of  the  magnetic  current  element  —  L1)-  However,  (8)  in  [9, 
sec.  22]  incorrectly  gives  the  previously  mentioned  /£  for  the  magnetic  field 
of  the  magnetic  current  element  nzS(z  —  l')-  This  error  is  pointed  out  in  [10] 
and  [11]. 

In  Appendix  D,  we  found  that  £(Q,i4<5(r  —  r')),  i4<5(r  —  r')),  j£(Q, 

y.z6(r  —  r')),  and  Jl(Q,ii26(r  —  r'))  are  given  by  (D.19),  (D.20),  (D.55),  and 
(D.58),  respectively.  In  these  equations,  k™ ,  7 ™ and  7^  are  given 
by  (B.7),  (B.24),  (B.41),  and  (B.53),  respectively. 
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Appendix  E 

Evaluation  of  Integrals  with 
Respect  to  4>'  and  <p 


In  Appendix  E,  the  integrals  (4.57)-(4.60)  and  (4.100)-(4.103)  are  evaluated. 


The  integrals  (4.57)-(4.60)  are 

*"  =  f‘%in(£^)Sin(r(#  -*'))<#'  (E.l) 

J<fo  0 

f2  =  [**  sin(^~— — )  cos(r(<£  -  <t>'))d<t>'  (E.2) 

J<t>  i  b 

<t>y3  =  [**  cos(^P)  sin(r(<£  -  <t>'))d<t>'  (E.3) 

J  <b\  b 

<fp4  =  f  cqs(~~—t — )  cos (r(4>  —  <f>'))d(f>'  (E.4) 

J<t>  l  o 


In  the  above  equations,  7  is  1  or  2,  y'1+  is  given  by  (4.39),  y'2+  by  (4.40).  <2>i 
by  (4.42),  and  <£2  by  (4.43). 

Changing  the  variable  of  integration  from  <f>'  to  y^+,  (E.1)-(E.4)  become 

4PX  =  /6sin(®^)sin(-(^+  -y7+)W+  (E.5) 

«/o  0  \X,j  / 

4>~'2  =  —  f  sin(^£ — )cos  (— (y'7+  -  y7+))  dy'7+  (E.6) 

x0  Jo  0  \X0  J 

<fp3  =  -- — —  f  cos(~ — )  sin  f  —  (y'7+  -  y7+)'j  dy'7+  (E.7) 

x0  Jo  b  \x0  ) 
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f4  =  —  f  cos(~——~  )  cos  f— (t/'y+  -  dyh 

X«  J 0  0  \X()  / 


(E.8) 


where  y7+  is  given  by  (2.15)  and  (2.16).  From  (1.2)  and  (2.18),  x0  is  given 

by 

*•-£  (E9) 
Using  [5,  formulas  401.02  and  401.04]  to  expand  the  trigonometric  functions 
of  the  difference  arguments  in  (E.5)-(E.8),  we  obtain 

-  ^lsin(^)}  (E.10) 

l  Xo  xo  J 


(jp2  =  <t>W  cos(  — — )  +  sin(— — ) 

X0  X0 

«r3  =  -(-ir(43>cos(^)  -  44>si 

l,  X0 

<p4  =  44)  cos(~ — )  +  43)  sin(— — ) 

Xo  x0 


■y  _i_  ^ 

sm(-— ) 
x0  . 


where 


4"  *  tJ  sin(?f)3iai7:)dx 

W  =  i  /\in(^)cos(— )dx 
p  x0  yo  o  x0 

<^,3)  =  —  f  cos(^y^)  sin(  —  )<fx 

X(j  «/o  O  Xq 

<^,4*  =  —  f  cos(~~~)  cos(— — ')<fx 

Xo  *<0  0  Xa 


(E.ll) 

(E.12) 


(E.13) 


(E.14) 
(E.15) 
(E.16) 
(E.  17) 


The  integrals  in  (E.14)-(E.17)  are  evaluated  by  using  the  integration  for¬ 
mulas  [5,  formulas  435,  445,  and  465] 

f  ■  ,  \  ■  ,L  1  f  sin((a  —  6)x)  sin((a  +  b)x)  | 

/  sin(ax)  sin(6x)dx  =  -  <  — - - -~J—L - — - ~—L  S  (E.18) 

J  2^0-6  a  +  6  J 

f  ,  ,  ,,  ,  ,  1  f sin((a  —  6)x)  sin((a  +  6)x)  1 

/  cos(ax) cos(6x) dx  =  -  { - -  -| - - - >  (E.19) 

J  2  {  a-  b  a  +  b  J 

[  .  ,  .  ..  .  ,  1  fcos((a-6)x)  cos((a  +  6)x)l  , _ 

/  sin(ax)  cos(6x)  dx  =  —  -  < - - — -  H - - - >  (E.20) 

9  2  (  a  —  o  a  +  o  J 


(E-1S) 


(E.19) 


(E.20) 
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If  |  a  |=|  b  |,  then  the  right-hand  sides  of  (E.18)-(E.20)  are  to  be  replaced 
by  their  limits  as  j  6 1  approaches  |  a|.  Unfortunately,  (E.20)  is  subject  to 
excessive  roundoff  error  when  either  (a  -  b)  or  (a  +  b )  is  small.  To  avoid  this 
error,  we  use  the  trigonometric  identity  [5,  formula  403.22] 

cos i  =  1  -  2sin2(|)  (E.21) 

to  recast  (E.20)  as 


J  sin(aa:)cos(6x)  dx  = 


2{^}  sin2 
'  + - 

a  —  b 


2  ^  (a+6)r  j 

a  +  b 


(E.22) 


The  “6”  which  appears  in  (E.18)-(E.20)  and  (E.22)  is  not  to  be  confused 
with  that  in  (E.14)-(E.17). 

Using  (E.18),  (E.19),  and  (E.22)  to  evaluate  the  integrals  in  (E.14)- 
(E.17),  we  obtain 


If  (pir  ±  ^)  is  zero,  then  the  right-Lacd  sides  of  (E.23)-(E.26)  are  to  be 
replaced  by  their  limits  as  (pir  i  ~-)  approaches  zero. 

The  integrals  (4.100)-(4.103)  are 


/  Vco. (23d w 

J  4>z  o 

,0,-72  (*'  .-72  •  /  rmrya+ 

4>  =  I  <jP  sm( - - - )d<p 

J  <b  3  0 
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•  ,Tn*ya+'J1 

<p  sin( - - — )a<p 

^3  0 

(E.29) 

4T*  =  / 

Jt 

'**  11 4  ,mxya+\ji 

<\T  cos( - - — )d<p 

h  0 

(E.30) 

where  7  is  1  or  2,  a  is  1  or  2,  y1+  is  given  by  (2.15),  y2+  by  (2.16),  <j> 3  by 
(4.104),  and  <£4  by  (4.105).  In  (2.15)  and  (2.16),  x0  is  given  by  (2.18)  where 
<f>0  is  given  by  (1.2).  In  (E.27)-(E.30),  <jP 1,  <p 2,  <pz,  and  <p 4  are  given  by 
(E.10)-(E.13),  respectively.  Substitution  of  (E.10)-(E.13)  into  (E.27)-(E.30) 
gives 


where 


irl  =  -(-1  r{4V272  -  4V271} 

(E.31) 

<r2  =  4V172  +  4V171 

(E.32) 

=  -(-1)7  {4V172  -  4V171} 

(E.33) 

=  <t>(f)<t>a2'r2  +  4V271 

(E.34) 

/ ai'yi  f*4  •  /miry°+x  •  /ry7+xji 

<f>  n  _  1  gm( - - )sm( - )<ty 

(E.35) 

b  X0 

(E.36) 

J<h  b  '  x0 

(E.37) 

±a  2-v2  /*4  ,  mirya+  v  ,rj^+ 

0  T  =  /  cos( - r — )cos( - )a0 

«/^3  0  X  ^ 

(E.38) 

Seeking  to  change  the  variable  of  integration  in  (E.35)-(E.38)  from  <p  to 
t/or+,  we  differentiate  (2.15)  and  (2,16)  to  obtain 


d<t>  =  ( 


dy 


Of  + 


dy 


X° 

a+ 


a  =  1 
a  =  2 


(E.39) 


Substituting  (4.104)  and  (4.105)  into  (2.15)  and  (2.16)  and  using  (E.9),  we 
obtain 


=  f  i>,  1  =  1 

J  4>=<t>3  (  0,  a  =  2 


(E.40) 
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»‘+U.  =  {“:“:2  <E-4» 

Applying  (E.39)-(E.41)  to  (E.35)-(E.38),  we  find  that 


X0  Jo  0 

■)sm(  )<fya+ 

x0 

(E.42) 

=  1  /‘cos(m7°+ 
x0  Jo  0 

-)sin(  )dt/a+ 

xQ 

(E.43) 

=  1  /‘sin  (m7“+ 
x0  Jo  0 

)cos(  )dya+ 

x0 

(E.44) 

>™  =  1  /*  cost  mTf+)  cost 

Xq  «/o  0  3/0 

a,  then  (E.42)-(E.45)  can  be  expressed  as 

(E.45) 

=  *£>, 

7  =  a 

(E.46) 

=  A31. 

7  =  a 

(E.47) 

7  =  a 

(E.48) 

v™  =  A1. 

7  =  c: 

(E.49) 

where  the  <^>’s  on  the  right-hand  sides  of  (E.46)-(E.49)  are  given  by  (E.14)- 
(E.17)  with  p  replaced  by  m.  If  7  ^  a,  we  add  (2.16)  to  (2.15)  in  order  to 
obtain 

j/7+  =  ttx0  +  b-  ya+,  7  /  a  (E.50) 

Substituting  (E.50)  into  (E.42)-(E.45)  and  using 
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we  obtain 


0171  =  (-l)r  < 

<f>m  Sin(~)  “  COS (~) 

{  xo  xo  J 

>  ,  7  #  cc 

(E.53) 

0271  =  (-1  )r  < 

4>^m  sin(~)  —  4$  cos  A 

L  x»  xo  J 

► ,  7  #  01 

(E.54) 

al 72  =  (-l)r  « 

<t>m  cos(““)  +  4$  sin(~) 

,  xo  xo  d 

S  7^a 

(E.55) 

0272  =  (-1)' . 

<t>m]  cosA  +  <j>£]  Bin  A 
[  xo  xo  J 

S  7  #  « 

(E.56) 

where  the  <f>' s  on  the  right-hand  sides  of  (E.53)-(E.56)  are  given  by  (E.14)- 
(E.17)  with  p  replaced  by  m.  These  <j>'s  can  be  calculated  from  (E.23)-(E.26). 

The  results  obtained  in  Appendix  E  are  stated  as  follows.  The  quantities 
(p1,  <£72,  jP3,  and  <jp4  are  given  by  (E.10)-(E.13)  in  which  <f>^\  and 

<j>W  are  given  by  (E.23)-(E.26).  The  quantities  <^an'1,  <£“72,  <£“73,  and  (j>ay4  are 
given  by  (E.31)-(E.34)  in  which  <f>£\  <j>(2\  4>f  \  and  <f>W  are  given  by  (E.23)- 
(E.26).  Moreover,  ,  <f>a2‘yl ,  2,  and  <f>a2‘y 2  are  given  by  (E.46)~(E.48) 

for  7  =  a  and  by  (E.53)-(E.56)  for  7^0. 
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Appendix  F 

Evaluation  of  Integrals  with 
Respect  to  J  and  z 


In  Appendix  F,  the  integrals  (4.77)-(4.88)  are  evaluated  for  (  —  |  <  z  <  |). 
Afterwards,  the  integrals  (4.106)-(4.110)  are  evaluated.  For  (-!  <  *  <  f)> 
the  integrals  (4.77)-(4.88)  are 


z“  =  I1  e-^-z')cos(SIllW  (F.l) 

J- !  c 

z62  =  J*  cosh  (7r,(^3  ~  *')  cos(^— )dz  (F.2) 

zTE3  =  jz  e-ilE{L3  -  z')  sin(i£f!t)<fc'  (F.3) 

zTEA  =  J*  sinh  (7-(X3  -  r'))  sin(^—  )dz  (F.4) 


where  z'+  is  given  by  (4.41).  In  (F.l)  and  (F.2),  8  is  either  T M  or  TE. 
Substituting 

7 6  =  TM,  TE  (F.5) 

in  (F.1)-(F.4),  using  [5,  formulas  654.6  and  654.7],  changing  the  variable  of 
integration  from  z'  to  z,+  ,  and  finally  replacing  z'+  by  x,  we  obtain 

Z*1  =  fo  {cos  (PSr,(L$  -  x)) -jsin  {$a{L$  -  x))}cos(^)dx  (F.6) 
ZS2  =  J ^  COS  ~  x))  cos (^— )dx  (F.7) 
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ZTEZ  -  Jo  {cos  (/3jE{L$  -  x))  -  j  sin  (/3rT4£(Lj  -  x))}  sin(^)dx  (F.8) 
zTEA  =  j  J+  sin  ({3j.E(Lt  —  x))  sin(^—  )dx  (F.9) 

where 


r+  =  z  + 1  (F.10) 

Lt=L  3  +  |  (F.ll) 


The  sum  of  (F.6)  and  (F.7)  is 

zsi  +  zS2  =  f  cos  —  x))  cos(- — )dx 

JO  c 

-j  sin(^4(IJ  -x))  cos(^)dx  (F.12) 


The  difference  between  (F.8)  and  (F.9)  is 

^TE3  _  ZT£4  _  J  CQS  ^TE(£+  _  sin^S-JU'jdx 

-j  jf  sin  -  x))  sin(^)dx  (F.13) 

Since  the  right-hand  side  of  (F.12)  is  simpler  than  that  of  (F.6),  we  will 
obtain  zsi  by  evaluating  ( z51  4-  zS2)  and  z52  and  by  setting 

z61  =  (zsi  +  zS2)  -  zS2  (F.14) 

Similarly,  we  will  obtain  zTE3  by  evaluating  ( zTE 3  —  zTE4)  and  zTE4  and  by 
setting 

zTE3  =  (zTE3  -  zTE4)  +  zTE4  (F.15) 

The  integrals  in  (F.7),  (F.12),  (F.9),  and  (F.13)  are  of  the  forms  (12. 
formulas  2.532] 


J  sin(ax  +  b)  sin(cx  +  d )  dx 


sin((a  —  c)x  4-  b  —  d) 

2  (a  —  c) 

sin((a  +  c)x  +  b  +  d) 
2  (a  4-  c) 


( F .16) 
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f  .  ,  .  ...  cos  ((a  —  c)x  +  b  -  d) 

/  sm  oi  4-  o) cosicz  -f  d)  dx  = - — - — 

J  2  (a  —  c) 


J  cos(aa:  +  6)  cos(cz  +  d)  dx  = 


cos((a  +  c)  x  +  b  +  d) 
2  (a  +  c) 

sin((a  —  c)x  +  b  —  d) 


(F.17) 


J  V  2(a_c) 

.  sin((a  +  c)x  +  6  +  d)  /r, 

+ - 2(a"+c) -  (F'18) 

In  (F.16)-(F.18),  a,  b ,  and  c  are  arbitrary  constants  not  to  be  confused  with 
the  specific  dimensions  a,  6,  and  c  in  Fig.  1.1.  If  a  =  ±c,  then  the  right-hand 
sides  of  (F.16)-(F.18)  are  to  be  replaced  by  their  limits  as  a  approaches  ±c. 

Using  (F.16)-(F.18)  to  evaluate  the  integrals  in  (F.7),  (F.12),  (F.9),  and 
(F.13),  we  obtain 

si  _  sin(g*-c  -f  0sraL$ )  -  sin (gs~z+  +  0sraL$) 

2  qs- 

sin(<75+c  -  0sraL$ )  -  sin {qs+z+  -  0srsL$ ) 

+ - 2^ -  (R19) 

si  .  52  sin (qs~c  +  PsrtL$)  ~  +  j  cos {qs~z+  +  0sraL$ ) 

2  +Z  =  - 2^~ - 

,  sin (qs+c  -  0srsL$ )  +  je~jPSr*Lt  -  j  cos {qs+z+  -  0sraL$ 

+ - -20) 


x*  = 


(F. 19) 


z01  +  z0i  = 


sin {qs+c  -  0sraL$ )  +  je  j^»Lt  -  j  cos(qs+z+  -  0sraL$ ) 


tea  -sin (qTE'c  +  QjE L$)  -  si n(qTE  z+  +  0jE L$ ) 

2  =_; - 2 7^ 

,  .sin (9T£+c  -  0jELt)  -  sin (qTE+z+  -  0™  L$) 

+J - -  (F-21) 

TE3  _ TEA  j  sin (qTE~c  +  PjfLt )  +  e~j@r3  Lt  -  cos (qTE~  z+  +  0jE  L$ ) 

2  "  2qTE~ 

-j  sin (gT£+c  -  0jEL$ )  +  e~iP™L$  -  cos {qTE+z+  -  0jE  L$ )  /T, 

- 2^r -  (F-22) 


(F.22) 


where 


q5~  =  31-0*  6  =  T  M,  T  E 


F.23) 
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<5  =  TM,TE 


( F -24) 


Substituting  (F.19)  and  (F  20)  into  (F.14)  and  using  (E.21),  we  arrive  at 


e  I  sin(q£~z+)  —  2j sin2(2^£^-) 

~~2  \  q*~ 

sin(q5+z+)  4-  2 j  sin2  (  9~y~  )  ] 

+  ?s+  } 


(F.25) 


Using  [5,  formulas  401.01  and  401.02]  to  expand  the  trigonometric  functions 
in  (F.19)  and  applying  (E.21),  we  obtain 


zS 2  =  [{  {sin(^  c)  -  sin(q*  z+)}  cos(0srsL$)  -  2{sin2(^y^) 
-sin2(^-^-)}sin(^,IJ)}/(295-)]  +  [{  {sin(q5+c)  -  sin(q5+z+)} 

•cos (ti.Lt)  +  2{sin2(^)  -  sin2(^)}  sin(/£L3+)}/(2q5+)]  (F.26) 

If  q8*  =  0,  then  the  right-hand  sides  of  (F.25)  and  (F.26)  are  to  be  replaced  by 
their  limits  as  q6±  approaches  zero.  If  0ST3  is  purely  real,  expressions  (F.25) 
and  (F.26)  are  suitable  for  calculating  (4.74),  (4.75),  (4.89),  and  (4.90)  at 
particular  values  of  z  because  (F.25)  and  (F.26)  are  not  subject  to  excessive 
roundoff  error  when  |q5±|  is  small. 

If  0sr3  is  purely  imaginary,  we  use  (F.5)  to  obtain 

=  -jit.  CF.27) 


Substituting  (F.27)  into  (F.25)  and  using  (E.21)  to  dispose  of  the  sin2  terms, 
we  arrive  at 


yS  /  + 


z&i  _  J± 


e~jqS  2+  _  i  l_ejq5+z+ ' 

-  +  - 7- -  > 


where 


qs  = 
?5+  = 


<F  .  s 
—  +  J7rJ 
c 

q*  ■  s 

--J  7r, 
c 


(F  .28) 


(F.29) 

(F.30) 
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Expression  (F.28)  is  recast  as 


:Sl  =  |e7r^+  ^ 


>.  ,qvz+,  q*  .  ,qirz+  ,N 

7r,  COs(- - )  +  ~  sin(~ - ) 

c  c  c  , 


-Is 

ITS 


e-frM 
(f)2  +  (7r5,)2 


(F.31) 


which  becomes 


zsi  =  j27*,e»7rj2  sinh(  )  +  e^raZ+ 


'qir  ._,q*z'¥s  s  .  2,qirz+ ^\]  < 


(F.32) 


Substituting  (F.27),  (F.29),  and  (F.30)  into  (F.19)  and  using  [5,  formula 
408.16],  we  find  that 


-Z52  =  "^7r»  cos(- - )  sinh  (7*,  (£3  -  *+))  -  —  sin(-^— —  )cosh 

c  c  c 

(7,S.(IJ  -  *+))  -  (-l)'7;,3inli  (-,‘r,(Ll  -  c))  }/  {(^)2  +  (%4,)2}  (F.33) 

If  both  ^  and  7^  are  zero,  then  the  right-hand  sides  of  (F.32)  and  (F.33)  are 
to  be  replaced  by  their  limits  as  7^  approaches  zero  while  q  is  held  at  zero. 
If  xTa  is  purely  real,  expressions  (F.32)  and  (F.33)  are  suitable  for  calculating 
(4.74),  (4.75),  (4.89),  and  (4.90)  at  particular  values  of  z  because  (F.32)  and 
(F.33)  are  not  subject  to  excessive  roundoff  error  when  q  is  zero  and  7^  is 
small. 

Substituting  (F.21)  and  (F.22)  into  (F.15)  and  using  (E.21),  we  arrive  at 
TE3  _  /  sip(<?rE~-?+)  ~  2J  si°2(  *TE2Z+  ) 

^  2  \  q™- 

sin (qTE+z+ )  +  2 j  sin2(  ) 

q™+ 

Using  [5,  formulas  401.01  and  401.02]  to  expand  the  trigonometric  functions 
in  (F.21)  and  applying  (E.21),  we  obtain 

zTE4  =  [{  {sin(?r£'c)  -  sin(/£_z+)} cos(^rr3£Z+) 


(F.34) 
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+2j {sin2( ^~2~ )  “  sin2(-  2  ~~)}  sin((3j3E L$)} /(2qTE~)\ 

+  [{;'  {sin(/£+c)  -  sin (9TE+^+)}  cos(^rrs£IJ ) 

+2;{sin2(^Y^)  -  sin2(  —  2  - ■ ) }  sin(/3jaEL$ ) }  /(2<?r£+)]  (F.35) 

If  =  0,  then  the  right-hand  sides  of  (F.34)  and  (F.35)  are  to  be  replaced 
by  their  limits  a s  qTE±  approaches  zero.  If  flJE  is  purely  real,  expressions 
(F.34)  and  (F.35)  are  suitable  for  calculating  (4.74),  (4.75),  (4.89),  and  (4.90) 
at  particular  values  of  z  because  (F.34)  and  F.oo)  axe  not  subject  to  excessive 
roundoff  error  when  |<frE±  |  is  small. 

If  0jE  is  purely  imaginary,  we  use  (F.5)  to  obtain 

*C  =  -hJ.E  (F.36) 


Substituting  (F.36)  into  (F.34)  and  using  (E.21)  to  dispose  of  the  sin2  terms, 
we  arrive  at 


ZTE3  _ 


e~7™Lt  fl 


•  TE-  + 
—  7/7  77 


-  e-jq 


„TE+r+ 


~TE- 


z'  l-ejq  z 

+  — 


,TE+ 


where 


qTB-=Sl+J1TB 

c 

TE+  _<]*__  -  TE 
H  J  irs 

C 

Expression  (F.37)  is  recast  as 

z™  =  {e^+  (7™,in (^)  -  ^cos(^)) 

+  c  J  (=)»  + 

which,  in  view  of  (E.21),  becomes 

=  {e^+  + 

\  c  c  2c  J 

-yTEL+ 

e  'rs  ^3 


2 qx  l~TB z+  .  7 JEz+^\ 

■  -t—ei'"*  smh (-^- — )  f 
c.  2  J 


(  =  )a  +  (7™)a 


(F.37) 

(F.38) 
( F .  39 ) 


(F  .40) 


(F.41 ) 
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Substituting  (F.36),  (F.38),  and  (F.39)  into  (F.21)  and  using  [5,  formula 
408.16],  we  find  that 


zTEA  =  {l™  sin(^-)  cosh(7rT/(Z+  -  2+)) 

+“  cos(^j~)  sinh(7 JaE(Lt  ~  z+)) 

-qA-iy  sinh (7rT/(Z+  -  c))}/{(^)2  +  (7rr/)2}  (F.42) 

c  c 

If  q  =  0,  then,  as  evident  from  (F.3)  and  (F.4),  the  right-hand  sides  of  (F.41) 
and  (F.42)  are  to  be  set  equal  to  zero,  regardless  of  the  value  of  7£E.  With 
this  reservation,  expressions  (F.41)  and  (F.42)  are  suitable  for  calculating 
(4.74),  (4.75),  (4.89),  and  (4.90)  when  7£E  is  purely  real. 

The  integrals  (4. 106)— (4. 1 10)  are 

*(”  =  f\  T™1  cosh(7™(£3  -  *)) 

+z™2e~1r*  ~  \  cos(  1~—Z'  ■  )dz  (F.43) 

J  c 

=  jf*  {2™  cosh(7;E(i3  -  *)) 

+:rBe- -  2> }  cos(— )dz  (F.44) 

J  c 

2(3)=  /*  {  _  sinh(7«(£3  _  ,)) 

J  ~~  2 

+zTB1t—tT,E{L3  -  3)}  sin(22l)<i2  (F.45) 

J  c 

*(i>  =  /_.  {^cosh(7-(i3--)) 

_2r£4e-7.r.E(L3  -  -)}  «*(  —  )<<.-  (F.46) 

J  C 

2<s>  =  /_*  i*™  -  ‘TE3>m-<lE(L3  - 

— rrE*e  ~  sin(  —  -  )d:  (F.47) 

J  c 
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where  z+  and  zTEs  are  given  by  (F.10)  and  (4.92),  respectively.  The  quanti¬ 
ties  z™1,  z™2,  zTEl ,  zTE 2,  zTE3,  and  zTE4  are  given  by  (F.1)-(F.4).  Chang¬ 
ing  the  variable  of  integration  from  z  to  z+,  substituting  j(3™  and  j/Sjf  for 
7™  and  ~iJE  in  accordance  with  (F.5),  and  using  [5,  formulas  654.6  and 
654.7],  we  recast  (F.43)-(F.47)  as 


i:  {(*“ 

+  z™7)  cos  0rM  -  jz™7  sin  0™  j  cos(  )dz 

J  C 

+(F.48) 

/; 

+  zTE2)  cos/?TE 

-  jzTE2 

sin  0TE 

\  <nrZ+  V A  + 

f  cos( - ) dz^ 

'  c 

(F.49) 

z<3>  = 

jf  {zTE2  cos  fE  -j(zTEl 

+  zTE2) 

sin  (3te 

1  .  ,n?rz+  + 
f  sin( - Jar 

J  c 

(F  .50) 

2W  = 

FM" 

TE4  \  qTE 

—  z  )  cos  p 

+  jzTE 4 

sin  fE 

\  <n1TZ+  + 
r  cos( - )dz 

J  c 

(F.51) 

r<5>  = 

FA*™- 

-  zTEt  coz  aTE  - 

j(zTE3- 

-  zTE 4) 

-'TC'I  .717 TZ 

sin/?r  r  sin( - 

J  c 

-)dr+ 

(F.52) 

where 

=  0i.(Li  -  ; 

O, 

6  =  TE 

,TM 

(F.53) 

in  which  Lj  is  given  by  (F.ll). 

Factors  appearing  in  (F.48)-(F.52)  are  [5,  formulas  401.05,  401.06,  and 
401.07] 


=  i{cos(n4+r+  -cos(n5_r+  +  #,££)}  (F.54) 

=  \{sin{ns+  z+  -  0sr,L$)  +  sin{ns~ z+  +  3sr3L$)}  (F.55) 
=  i{  sin(n5"2+  +  $tL$)  -  s\n{ni+ z+  -  (F.56) 

=  cos{ns~z+  +  +  cos{n6+z+  -  (F.57) 

where 


mrz 


sin  /?4  sin( - ) 


mrz 


cos  (3 4  sin( - ) 


mrz 


sin/?4  cos( - ) 


mrz 


cos  /?4  cos( - ) 


n4~  = 


ns+  = 


-  -  #, 

c 

TllT  r,t 

—  +  A?r 

c 
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( F  .58) 
( F .59) 


in  which  8  is  either  TE  or  TM.  Substituting  (F.19),  (F.20),  (F.56),  and 
(F.57)  into  (F.48)  and  (F.49),  we  obtain 

z(1)  =  z™  (F.60) 

z{2)  =  zTE  (F.61) 

where 

=  [{(sin(?‘-c  +  tfaLt)  ~  )(zsc. _  +  zse+)  -  j(zsa_  -  zsa+) 

•sin  {qs~c  + pSraL+)  +  >(*?_,,_  +  zse_ie_)  -  j(zsa+a_  -  Ze+,c-)}/(4g4-)] 

+  [{(sin(?5+c  -  psraL$ )  +  je~iPr,Lt  )(2f_  +  Z6C+)  -  j(zsa_  -  z6a+) 

•sin {qS+c  -  ~  2*-,c+)  “  i(zU,*+  +  4+,e+)}/(V+)]  (F-62) 

in  which 

zaa  =  I  sin(nfaz+  -  (al)tf,L$)dz+  (F.63) 

Jo 

*L=  f  cos(n‘"2+  -  (a\)P‘„LS)dz+  (F.64) 

JO 

*L„  =  -  (al)tf,i;)Sin(q‘V  -  (7I)rf,l+)d;+  (F  65) 

zL,cy  =  I  cos{n5az+  -  {al)l3sraL$)cos(qHz+  -  {-fl)j3sraL$)dz+  (F.66) 
Jo 

Here,  al  =  +1  when  a  =  +,  and  al  =  —1  when  a  =  — .  The  quantity  7I 
is  similarly  defined.  Substitution  of  (F.19),  (F.20),  (F.54),  and  (F.55)  into 
(F.50)  yields 

z<3>  =  [{(4*  +  &«*?*-*  +  **«)  -  U  »(«”-«  +  JlELi ) 
+e-jt?,ELl  )(ZTB  _  ££)  _  (ZTE'_  _  zTEc  _  (£B  _  +  .T*.)}/^™-)] 

+  [{<*"  +  *,T-£)  sia(?T£+C  -  0lEL+)  +  (-;  sm(,rE+c  -  0™ L$ ) 

+e-jAr.£4 )(££  _  2T£)  _  +)  -  (z™,  +  -  .-fE .,)}/( V£+)] 

(F.G7) 
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Evaluating  the  integrals  in  (F.63)  and  (F.64),  we  obtain 


ZL  = 


cos {P6rtL$)  -  cos {nSac  -  (al )P*,L$) 


s  sin (nSac  -  (al )/3sraL$ )  +  (al)sin {ftM ) 


** COt 


(F.68) 

(F.69) 


If  nSa  is  zero,  then  the  right-hand  3ides  of  (F.68)  and  (F.69)  are  to  be  replaced 
by  their  limits  as  nSa  approaches  zero.  Application  of  the  integration  formula 
(F.16)  to  (F.65)  gives 


s  sin((n5a  -  g^)c+  (7I  -  al -  sin((7l  -  al)(3sr3L$) 
Z'a’r'  2{nSa  -  q*) 

sin((n5a  +  q*)c  -  (al  +  7 +  sin(M  + 

2  (nSa  +  qs~') 

Application  of  the  integration  formula  (F.18)  to  (F.66)  gives 


(F.70) 


- 


sin((nga  -  gs~')c  +  (7I  -  al)/3sraL$)  -  sin((7l  -  al)0fsL$) 

2 (nSa  -q*>) 

sin((nga  gs^)c  —  (al  +7 l)/3*aL$)  +  sin((orl  + -fl)(36r,LZ) 

2{nSa  +  q*) 


(F.71) 


The  sum  of  (F.70)  and  (F.71)  is 


&  sin ((nSa  -  <?*7)c  4-  (7I  -  al)ffltL$)  -sin((7l  -  al)$tL$) 


z  +  z 

3a, ry  1  ca,cry 


(n6a  -  q *) 


( F  .72) 


The  difference  between  (F.70)  and  (F.71)  is 


s  _  s  sin{{nSa  +  gs',)c  -  (al  Hr  7I ±  sin((al  +  ~/l)35rsL^ 

z,a,ry  Zca,<n  (n5a-|-^) 


(F.73) 
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Substitution  of  (F.68),  (F.69),  (F.72),  and  (F.73)  into  (F.62)  and  (F.67) 
gives 


=  ie-j^rsLt 
4 

r  (e~Jw*  c  -l)sm(qs~c  +  f3^Lt)  -  sin(n*~c  +  ftaL$)  +  sin(/36T3Lt) 

(1  ~ejn  +c)sin(g5~c  +  PsrsL%)  -  sin (ns+c  -  @*aL$)  -  sin 

ns+qs~ 

(e-jn  c  _  L)  sj n(qs+C  -  ffaL$)  +  sin (ns~c  +  PrM)_  ~ 

ns~qs+ 

(1  ~  einS+c)sm{gs+c  -  +  sin (n*+c  -  +  sin 

ni+g*+  / 

j  f  sin((n^~  -  gs~)c)  sin((n*+  +  <?*~)c) 

4  *■  (ns~  —  (n5+  +  qs~)qs~ 

sin((ns“  +  (f+)c)  sin ((ns+  —  \ 

(ns~  +  qs+)qs+  (ns+  -  qs+)qs+  * 
j(3)  .  \e-J^Lt 


(F.74) 


(1  -  e  Jn  c)  sin(gT£  c  +  &*EL% )  +  sin(nT£  c  +  0*ELt )  -  sin {/3jE  ) 


nTE-qTE- 


(1  -eJ»  +c)  sin(qrE~c  +  f3jE L+ )  _  sin(nr£+c  -  ~  sin(ffg/;+) 


nTE+qTE- 


(1  -  e  7n  c)  sin(<?r£+c  -  /3rT/l^ )  ~  sin(nTE  c  +  #T£ I3  )  +  sin^/l^ ) 


nTE-qTE+ 


(1  -  e-?n  c)  sin((7T£+c  -  )  +  sin (nTE+c  -  (3jEL$ )  +  sin [BjEL$ ) 


nTE+qTE+ 


If  sin ((nTE~  —  qrE  )c)  sin((nr£+  +  qTE  )c) 
4  ^  (nTE~  —  qTE~)qTE~  (nTE+ +  qTE~)qTE~ 


sin  ((nTE  +qTE+)c)  sin  ((nTE+ —  qTE+)c) 
'  (nTE-  +  qTE+}qTE+  ~  (nT£+  _  qTE+)qTE+ 


( F .75) 
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The  term  {sin((n5~  —  9i-)c)}/{(n'5-  —  qs~)qs~}  in  (F.74)  becomes  infinite 
as  qs~  approaches  zero.  We  subtract  {sin(n5-c)}/(n5-^-)  from  this  term  to 
render  it  finite  as  qs~  approaches  zero.  To  maintain  equality  in  (F.74),  we  add 
an  appropriate  quantity  to  the  term  whose  denominator  is  ns~q5~.  Similarly 
treating  the  ill-behaved  terms  whose  denominators  are  (n5+  +  qs~)q°~ ,  ( n 8  + 
<75+)?i+i  and  (n5+  —  (^+)qs+  in  (F.74),  we  arrive,  with  the  help  of  [5,  formulas 
401.01  and  401.02],  at 

zs  =  J-(D*GS  +  c2Fs)  (F.76) 


where 


in  which 


D6  = 
G8  = 

F*  = 


1  ,  .-iBLlA 


n 


s- 


n 


6+ 


sin(<75~c  +  )  ~  sin {j%M_ ) 

q8' 

.  ain (qs+c  -  06raL$ )  +  sm{f3srtL$ ) 

+  ^ 

-f(n6~c,  - qS~c )  +  /(n5+c,  qs~c) 

-/(n5~c,  qS+c)  +  /(n*+c,  -qs+c) 


f{x 


sin(x  4-  y)  sinx 


x  +  y 


(F.77) 


(F.78) 

(F.79) 

(F.80) 


If  we  divided  the  right-hand  side  of  (F.74)  by  j,  changed  the  signs  of  the 
terms  whose  denominators  are  n5~qs~,  ( n6~  —  qs~)q&~ ,  n6~qs+ ,  and  (n5-  + 
qs+)qs+ ,  and  then  replaced  6  by  TE,  we  would  have  the  right-hand  side  of 
(F.75).  Therefore,  we  can,  from  inspection  of  (F.76)-(F.80),  write 


x<3>  =  i(£(3)G(3)  +  c2F<3)) 


(F.81) 


where 


D{3) 

g(3) 


=  G 


TE 


(F.S2) 

(F.83) 
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(F.84) 


F(3)  =  f(nTE  c,-qTE~c)  +  f(nTE+c,qTE  c ) 

+f(nTE~c,  qTB+c )  +  f(nTE+c,  - qTE+c ) 

In  (F.83),  GTE  is  given  by  (F.78)  with  6  replaced  by  TE. 

In  this  paragraph  and  the  next  two  paragraphs,  we  obtain  formulas  that 
are,  when  is  purely  real,  suitable  for  calculating  D s,  D^3\  Gs ,  and  the  /’ s 
in  (F.79)  and  (F.84).  Concerned  about  roundoff  error  when  |n5±|  is  small, 
we  use  (E.21)  to  express  D6  of  (F.77)  and  of  (F.82)  a s 


Ds  =  { 

£><3>={ 


,ns~c 


n5+c. 


— jsin(ns~c)  —  2  sin2 (  —  — )  —  j  sin(n's+c)  +  2  sin2(— 7 — ) 


&  r  + 


n 


s- 


+ 


n 


6+ 


(F.85) 


,  nTE~c 


nTE+c . 


j  sin(nr'B-c)  +  2sin2( — - — )  —  j  sin(nT£’+c)  +  2sin2( — - — ) 


n 


TE- 


n 


TE+ 


} 


(F.86). 


If  ns±  —  0,  then  the  l/ns±  terms  in  (F.85)  and  (F.86)  must  be  replaced  by 
their  limits  as  ns±  approaches  zero.  The  right-hand  sides  of  (F.85)  and  (F.86) 
were  purposely  expressed  so  that  the  values  of  these  limits  are  obvious. 

Concerned  about  roundoff  error  when  \qs±  |  is  small,  we  use  [5,  formulas 
401.01  and  401.02]  and  (E.21)  to  express  Gs  of  (F.78)  as 

sin {qs~c)  cos -  2  sin2  {— -)  sin (#,/£) 

G  - - - - 2 - 

sin (qs+c)  cos(^s£^ )  +  2  sin2  (^—)sin (0*raL$) 

+ - ^ - 2 -  (F.S7) 

If  qs±  =  0,  then  the  right-hand  side  of  (F.87)  must  be  replaced  by  its  limit 
as  q5±  approaches  zero. 

Caxe  must  also  be  taken  to  avoid  excessive  roundoff  error  in  the  calcula¬ 
tion  of  f{x,y )  of  (F.80).  From  (F.23),  (F.24),  (F.58),  and  (F.59),  we  have 

(n5_  -  qs~)c  =  (n  -  9)tt 
(ns+  +  q6~)c  =  (n  +  q)  n 
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(F.88) 

(F.S9) 


(n5  +  qs+)c  =  (n  +  q)n  (F.90) 

(n5+  —  qS+)c  =  (n  —  q)n  (F.91) 

so  that,  for  every  /  that  appears  in  (F.79)  and  (F.84), 

x  +  y  =  (n  ±  q)ir  (F.92) 

Equation  (F.92)  reduces  (F.80)  to 

sin  t 

f(x,y)  = - ,  x  +  y^Q  (F.93) 

yx 

f(x,y)  =  - — 7—^,  x  +  y  =  0  (F.94) 

y 


Expression  (F.94)  was  obtained  by  taking  the  limit  of  the  right-hand  side  of 
(F.80)  as  x  +  y  approaches  zero.  If  |y|  <  f ,  we  use  (F.92)  to  recast  (F.93)  as 


f{x,y) 


(_l)n*q  siny  f  i  +  y/  0 
yx  ’  1  |y|  <  f 


(F.95) 


If  |yj  <  0.1,  we  replace  the  right-hand  side  of  (F.94)  by  the  series  approxi¬ 
mation  [5,  formula  415.01] 


y  v3  v5 

/(*’»)  =  31  -  51  +  v.  ' 


x  +  y  =  0 

|y|  <0.1 


Collecting  the  results  (F.93)-(F.96),  we  have 


f{x,y)  =  { 


sin  x 


_  j  x  +  y  #  0 

yx  \  |y|  >  f 

(~l)n±l?  sin  y  (  x  +  y  ±  0 
yx  \  |y|  <  f 

y  —  sin  y 


(  x  +  y  =  0 

\  |y|  >  0-1 

l _yl  ,yl  I  x+y  =  0 

.3!  5!  +  7!’  \  M  <0.1 


(F.96) 


(F.97) 


where  (n  ±  q)  is  the  integer  that  satisfies  (F.92). 
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In  this  paragraph  and  the  next  paragraph,  we  obtain  formulas  that  are, 
when  (36ra  is  purely  imaginary,  suitable  for  calculating  Ds,  Gs ,  Fs,  and 

F^3\  We  arrive  at  these  formulas  by  substituting  (F.27),  (F.29),  (F.30),  and 


ns~  = 


n‘+  = 


nir  .  tt 

—  +  nL 

c 

™  .  s 

~  ~ ]lra 


(F.98) 

(F.99) 


into  (F.77),  (F.82),  (F.78),  (F.79),  and  (F.84).  Substitution  of  (F.27),  (F.98), 
and  (F.99)  into  (F.77)  gives 


D6 

Ds 


*^sinh(-^^) 
(titt)2  +  (7  srac)2 

4j cosh(-^p) 
(nit)2  -I-  (7rsc)2 


n  even 


n  odd 


(F.100) 

(F.101) 


If  both  n  and  7*,  are  zero,  then  the  right-hand  side  of  (F.100)  must  be 
replaced  by  its  limit  as  7^  approaches  zero  while  n  is  held  at  zero.  This  limit 
is  —2 jc.  Substitution  of  (F.27),  (F.98),  and  (F.99)  into  (F.82)  gives 


TE 

Anxc  -f)  sinh(^-~) 

n(3)  _ _ 2 

(nit)2  +  (~fJ,Ec)2 

TE 

4n?rc  cosh(  — - ) 

n(3)  _  _ 2 _ 

(tit)2  +  (7rr,Ec)2 


n  even 


n  odd 


(F.102) 

(F.103) 


If  both  n  and  are  zero,  then  the  right-hand  side  of  (F.102)  must  be 
replaced  by  its  limit  as  approaches  zero  while  n  is  held  at  zero.  This 
limit  is  zero.  Substitution  of  (F.27),  (F.29),  and  (F.30)  into  (F.78)  and  use 
of  [5,  formula  408.16]  lead  to 


s  _  27ry  {sinh(7r*al;t)  -  (— l)*sinh(7 -  c))} 
(qit)2  -I-  (7^,c)2 


(F.104) 


If  both  q  and  7®,  are  zero,  then  the  right-hand  side  of  (F.104)  must  be  replaced 
by  its  limit  as  7*,  approaches  zero  while  q  is  held  at  zero.  This  limit  is  2c. 
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Substituting  (F.80),  (F.29),  (F.30),  (F.98),  and  (F.99)  into  (F.79),  using 
[5,  formula  408.16],  and  assuming  that  n  is  precisely  an  integer,  we  obtain 

pS  -  fsin((n-j-g)x)  sin((n  -  q) x) 

(qn)*  +  (7r,c)2  (n  +  q)*  (n  ~  q)* 

27«c(-l)n  sinh(7r54c)  \  /1?  in=x 

■  («^  +  R.c)>  /  (F'105) 


Similarly,  (F.84)  yields 

r,(3) _ 2x  f  /sin((n  +  g)x)  _  sin((n  -  g)xK 

(q*)2  +  ('lTtEc)2'*q'  ( n  +  q)ir  (n  -  q)ir 

2n7rT/g(-l)n  sinh(7rr/c)  \ 
(nx)2  +  (7J/C)2  J 


(F.106) 


If  q  —  ±n,  then  the  right-hand  sides  of  (F.105)  and  (F.106)  must  be  replaced 
by  their  limits  as  q  approaches  ±n  while  n  is  held  at  its  integer  value.  Thus, 
assuming  that  both  q  and  n  are  non-negative  integers,  forms  suitable  for 
calculation  are: 


f  q  =  n  =  0 

\  7r»C  >  0.1 


47  ( 7rf,c)2(-  l)n  sinh(7rssc) 

{W2  +  (7rV)2}{(nx)2  +  (7^c)2}’. 

2j7r»c  f  27r»c(— l)n  sinh(7r4c) 

M2  +  (7raC)2  1  M2  +  (7r,c)2 
4;  {7r,c  -  sinh(7r5,c)}  [  q  =  n  =  0 
(7rac)2  ’  l  7r6,C  >  0.1 

f7r,C  .  (7r*C)3  ,  (7rsC)5  I  f  <Z  = 

4Tir  +  -5r  +  -^rr  w,c 

47rn7r»E  c(~  l)n  sinh(7rT4Ec) 

{(gx)2  +  (7™c)2}  {(nx)2  +  (7rr/c)2}  ’ 

2xn  f27^Ec(-l)nsinh(7j/c) 

(9x)2  +  (7 JaEc)2  \  (nx)2  +  (7 JsEc)2 

0,  q  =  n  =  0 


q  =  n  ^  0 


■  f  7r,c  ,  (7*,<03  (t£c)«1  f  ?  =  n  =  0 
?  \  3!  +  5!  +  7!  J  ’  \  7r,c  <  0-1 


9  ±  n 


(nx)2  +  (7rr/c)2 


l|  ,  q  =  n 


(F.107) 

(F.108) 

(F.109) 

(F.110) 

(F.lll) 

(F.112) 

(F.113) 


If  q  =  0  or  n  =  0  and  if  7 ^  =  0,  then  (F.107)  and  (F.lll)  must  be  replaced 
by  their  limits  as  7*,  approaches  zero.  The  series  expansion  [5,  formula  657.1] 
was  used  to  obtain  (F.110). 
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(F.114) 


Comparing  (F.21)  with  (F.19),  we  note  that 

gTE4  _  •  -TE2 

where  zTE2  is  the  right-hand  side  of  (F.21)  with  the  sign  of  the  1  /qs+  term 
changed  and  with  6  replaced  by  TE.  Comparing  (F.22)  with  (F.20),  we  find 
that 

ZTEZ  _  ZTE4  =  j^TEl  +  5TE2)  (F .115) 

where  ( zTEl  +zTE 2)  is  the  right-hand  side  of  (F.20)  with  the  sign  of  the  l/qs+ 
term  changed  and  with  6  replaced  by  TE.  Substituting  (F.114)  and  (F.115) 
into  (F.51)  and  (F.52)  and  comparing  with  (F.49)  and  (F.50),  we  obtain 

(F.116) 

(F.117) 

where  z^  is  the  right-hand  side  of  (F.49)  with  the  signs  of  the  l/qTE+  terms 
changed.  Similarly,  i*3)  is  the  right-hand  side  of  (F.50)  with  the  signs  of  the 
lfqrE+  terms  changed.  According  to  (F.61),  z W  is  given  by  the  right-hand 
side  of  (F.76)  with  S  replaced  by  TE. 

Substituting  (F.76)  into  (F.116)  and  noting  that  the  integrations  that 
were  performed  in  obtaining  (F.76)  from  (F.49)  did  not  introduce  any  1  /qs+ 
factors  in  addition  to  those  in  (F.19)  and  (F.20),  we  obtain 

z(4)  =  -l(DTEG{4)  +  c2F(4))  (F.118) 

where  DTE  is  given  by  (F.77)  with  6  replaced  by  TE,  G(4)  is  given  by  the 
right-hand  side  of  (F.78)  with  the  sign  of  the  1  /qs+  term  changed  and  with  6 
replaced  by  TE,  and  is  given  by  the  right-hand  side  of  (F.79)  with  the 
signs  of  the  1  /q6+  terms  changed  and  with  S  replaced  by  TE.  Substituting 
(4.92)  and  (F.81)  into  (F.117),  noting  that  the  integrations  that  were  per¬ 
formed  in  obtaining  (F.81)  from  (F.50)  did  not  introduce  any  l/qTE+  factors 
in  addition  to  those  in  (F.19)  and  (F.20),  and  using  (F.5),  we  obtain 

. nTE  TE 

*(5)  =  fuTkii  +  +  c2F(5))  (F.119) 

{“ra  )  ** 


*<*>  =  jJ2 

z™  ,m(—)dz+ 
Jo  c 
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where  is  given  by  (F.82),  is  the  same  as  in  (F.118),  and  F <5)  is  given 
by  the  right-hand  side  of  (F.84)  with  the  signs  of  the  l/qTE+  terms  changed. 
Moreover, 


z™  =  /  sin(^-^— )sin (^^—)dz+  (F.120) 

Jo  c  c 

Assuming  that  both  n  and  q  are  non-negative  integers,  we  have 


f  f,  n  =  ?#0 

\  0,  otherwise 


(F.121) 


If  0^E  is  purely  real,  then  DTE  and  are  suitably  given  by  (F.85) 
and  (F.86).  Assuming  that  0*E  is  purely  real,  we  proceed  to  obtain  suitable 
expressions  for  G^4\  F^4\  and  F^s\  From  (F.87),  we  have 


<?E~c 


sm{qrE-c)cos{0jaELt)-2sin2{ 

G{4)  - - 2 


)  sin  {0*EL$ ) 


~TE— 


sin(qrE+c)  cos  (0™  L$ )  +  2  sin2(  —  - -C)  sin  (0?EL$ ) 

— - 


(F.122) 


From  (F.79)  and  (F.84),  we  have 

F{4)  =  -  f{nTE'c ,  -qTE-c)  +  f(nTE+c,  qTE~c ) 

+f{nTE~c,  qTE+c)  -  f{nTE+c ,  -qTE+c)  (F.123) 
F<5>  =  f(nTE-c,  - qTE~c )  +  f(nTE+c,  qTE~c) 

-f(nTE'c,  qTE+c )  -  f{nTE+c ,  - qTE+c )  (F.124) 


where  f{x,y)  is  given  by  (F.97). 

If  (3?E  is  purely  imaginary,  we  point  out  that  DTE  is  suitably  given  by 
(F.100)  and  (F.101)  with  6  replaced  by  TE  and  that  D ^  is  suitably  given 
by  (F.102)  and  (F.103).  Assuming  that  /3jE  is  purely  imaginary,  we  proceed 
to  obtain  a  suitable  expression  for  G^.  From  (F.78),  we  have 


G(4)  _  sin (gTE'c  +  f3jEL$ )  -  s\n((3jaE L $ ) 

?TE_ 

sin(qTE+c  -  0jEL%)  +  sin (0jtEL$) 
qTE+ 


(F.125) 
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(F.126) 


Substitution  of  (F.27),  (F.29),  and  (F.30)  into  (F.125)  leads  to 

(4)  _  j2Trgc{sinh(~tTaEL$)  -  (-1)*  sinh(7%g(Z£  -  c))} 

(qir)2  +  ilJaEc)2 

If  both  q  and  ~ftE  are  zero,  then  the  right-hand  side  of  (F.126)  must  be 
replaced  by  its  limit  as  7 JE  approaches  zero  while  q  is  held  at  zero.  This 
limit  is  zero. 

In  this  paragraph,  formulas  are  obtained  for  calculating  Fty  and 
when  (iJE  is  purely  imaginary.  Substituting  (F.29),  (F.30),  (F.98),  and  (F.99) 
into  (F.123),  using  (5,  formula  408.16],  and  assuming  that  n  is  precisely  an 
integer,  we  obtain 


p(4)  = 


2r  q  f  sin((n  +  q)i r)  sin((n  —  q)ir) 

TEF 2  L  (n  +  q)ir 


(qr)2  +  (~tJEc) 


(n  +  q)ir  '  (n  -  q)x 
27rr/c(~^)n  sinh(7^/c) 
(™02  +  {iJEc)2 


} 


(F.127) 


Similarly,  (F.124)  yields 


F<5)  = 


2  j 


te  /»in((n  +  ?)y)  sin((n  -  q)*)  \ 

{qir)2  +  (7 JEc)2  ^  7r*  '  (n  +  q)n  (n  -  q)v 

,  2n2nq(—l)n  sinh(7 JEc)  j  ^ 


(htt)2  +  (7^/c)2 


If  q  =  ±n,  then  the  right-hand  sides  of  (F.127)  and  (F.128)  must  be  replaced 
by  their  limits  as  q  approaches  ±n  while  n  is  held  at  its  integer  value.  Thus, 
assuming  that  both  q  and  n  are  non-negative  integers,  forms  suitable  for 
calculation  are: 


p(A)  =  47r<Z7rr/c(-l)n  sinh(7rTa£c)  , 

{{qir)2  +  (7 JsEc)2)  {(n7r)2  +  (7rT,Ec)2}  ’ 
p(4)  2 vq  f  _  27rT/c(-l)n  sinh(7rT5£c)  \ 

(q*)2  +  {'iJtEc)2  ^  (mr)2  +  (7 JEc)2  * 

F(4)  =  0,  q  =  n  =  0 
p( 5)  _  4jffang(- l)n  sinh(7rT/c) 

M2  +  (7 ’  q  * 


(F.129) 


q  =  n  ^  0  (F.130) 
(F.131) 
(F.132) 
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F(S)  2j  fTE  2(n7T)2(-l)n3inh(7rr/c)\ 

(?*)2  +  (7 lEcY  Xlrs  +  (nx)»  +  (7«c)’  J 

=  0,  <7  =  n  =  0 


(F.133) 

(F.134) 


Ifg  =  0orn  =  0  and  if  =  0,  then  the  right-hand  side  of  (F.129)  must 
be  replaced  by  its  limit  as  'yjf  approaches  zero.  If  n  =  0  and  j^tE  =  0,  then 
the  right-hand  side  of  (F.132)  must  be  replaced  by  zeroi. 

The  results  obtained  in  Appendix  F  are  cataloged  in  Table  F.l.  In  Ta¬ 
ble  F.l,  the  quantity  in  the  first  column  is  given  by  the  equatiou  whose 
number  appeaxs  in  the  second  column  when  the  nature  of  either  1™  or  7rT/ 
is  indicated  in  the  third  column. 
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Table  F.l:  Results  in  Appendix  F 


Quantity 

Equation  Number 

Circumstance 

mam 

(F.25) 

muni 

(F.26) 

purely  imaginary 

BliBH 

(F.32) 

(F.33) 

Kiapmamtmm 

EH 

(F.25) 

purely  imaginary 

BiEBP 

(F.26) 

B 

(F.34) 

7^b'  purely  imaginary 

(F  .35) 

I™  purely  imaginary 

(F.32) 

purely  real 

(F.33) 

'yjf  purely  real 

(F.41) 

-^TE4 

(F.42) 

7^E  purely  real 

EH 

(F.60),(F.76),(F.85),(F.87), 

(F.79),(F.97) 

H&lf.rHBIHnHHBHnMHH 

EU|m 

&■ 

(F.60),(F.76),(F.100),(F.101), 

(F.104),(F.107)-(F.110) 

7™  purely  real 

mm 

(F.60),(F.76),(F.85),(F.87), 

(F.79),(F.97) 

7^E  purely  imaginary 

mm 

(F.81),(F.86),(F.83), 

(F.87),(F.84),(F.97) 

7 JSE  purely  imaginary 

HB 

(F.118),(F.85),(F.122), 

(F.123),(F.97) 

7^b  purely  imaginary 

mmn 

(F.119),(F.121),(F.86), 

(F.122),(F.124),(F.97) 

7^E  purely  imaginary 

tm 

(F.60),(F.76),(F.100),(F.101), 

(F.104),(F.107)-(F.110) 

7^E  purely  real 

wm 

(F.81),(F.102),(F.103),(F.83) 

(F.104),(F.111HF.113) 

7EE  purely  real 

mm 

(F.118),(F.100),(F.101), 

(F.126),(F.129)-(F.131) 

7E£  purely  real 

VM 

(F.119),(F.121),(F.102), 

(F.103),(F.126),(F.132)-(F.134) 

7EE  purely  real 
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